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LATTICE BASED MULTISPACES 


GLORIA MAR{ BEFFA AND ELIZABETH L. MANSFIELD 


ABSTRACT. In this paper, we develop the theory of the discrete moving frame 
in two different ways. In the first half of the paper, we consider a discrete 
moving frame defined on a lattice variety, and the equivalence classes of global 
syzygies that result from the first fundamental group of the variety. In the sec- 
ond half, we consider the continuum limit of discrete moving frames as a local 
lattice coalesces to a point. To achieve a well-defined limit of discrete frames, 
we construct multispace, a generalization of the jet bundle that also gener- 
alizes Olver’s one dimensional construction. Using interpolation to provide 
coordinates, we prove that it is a manifold containing the usual jet bundle as 
a submanifold and show that continuity of a multispace moving frame ensures 
that the discrete moving frame converges to the continuous one as lattices co- 
alesce. We prove further that the discrete invariants and syzygies approximate 
their smooth counterparts. In our last chapter we discuss two important ap- 
plications, one to the discrete variational calculus, and the second to discrete 
integrable systems. Finally, in an appendix, we discuss a more general result 
concerning equicontinuous families of discretisations of moving frames, which 
are consistent with a smooth frame. 
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1. INTRODUCTION 


The theory and the applications of Lie group based moving frames is now well 
established, and provides an “invariant calculus” to study differential systems which 
are either invariant or equivariant under the action of a Lie group. Associated with 
the name of Elie Cartan [10], who used repéres mobile to solve equivalence problems 
in differential geometry, the ideas go back to earlier works, for example by Cotton 
[13] and Darboux [14]. 

Moving frames were further developed and applied in a substantial body of 
work, in particular to differential geometry and exterior differential systems; see 
for example papers by Green [19] and Griffiths [20]. From the point of view of 
symbolic computation, a breakthrough in the understanding of Cartan’s methods 
for differential systems came in a series of papers by Fels and Olver [16, 17], Olver 
(50, 51], Hubert [27, 28, 29], and Hubert and Kogan [30, 31], which provide a 
coherent, rigorous and constructive moving frame method. The resulting differential 
invariant calculus is the subject of the textbook, [37], and there have been an 
extensive number of applications since its publication, including to the integration 
of Lie group invariant differential equations, to the Calculus of Variations and 
Noether’s Theorem, (see for example [36, 21, 22]), and to integrable systems (for 
example [39, 42, 43, 44]). 

The first results for the computation of discrete invariants using group-based 
moving frames were given by Olver who called them joint invariants in [51]; mod- 
ern applications to date include computer vision [49] and numerical schemes for 
systems with a Lie symmetry [12, 33, 34, 35, 40]. While moving frames for discrete 
applications as formulated by Olver do give generating sets of discrete invariants, 
the recursion formulae for differential invariants which were so successful for the 
application of moving frames to calculus based results, do not generalize well to 
joint invariants. In particular, joint invariants do not seem to have recursion for- 
mulae under the shift operator that are computationally useful. To overcome this 
problem, the authors, together with Jing Ping Wang, introduced the notion of a 
discrete moving frame which is essentially a sequence of frames [41]. In that pa- 
per we prove discrete recursion formulae for a small computable generating sets of 
invariants, which we call the discrete Maurer—Cartan invariants, and investigated 
their syzygies, that is, their recursion relations. The main application to date has 
been to discrete integrable systems, with the authors of [46] proving that discrete 
Hamiltonian structures for W,,-algebras can be obtained via a reduction process. 
We note that a sequence of moving frames was also used in [34] to minimize the 
accumulation of errors in an invariant numerical method. 

In this paper, we extend the theory of discrete moving frames in two ways. The 
first is to consider a discrete moving frame defined on a lattice variety, which can be 
thought of as the vertices, or 0-cells, together with their adjacency information, in 
a discrete approximation of a manifold. We describe their associated cross-sections 
and define Maurer-Cartan invariants and local syzygies. In Section 3.2 we further 
classify global syzygies and prove that they are associated to topological aspects 
of the variety, like representatives of the discrete fundamental group of the lattice 
variety, with properties like twisting. 

The second extension, beginning in Section 4 and for the rest of the paper, is 
to consider families of discrete frames and how their continuum limits may define 


DISCRETE MOVING FRAMES 3 


smooth frames. Our interest in this second case is how discrete invariants and their 
recursion relations limit to differential invariants and their differential syzygies. 

In order to provide a general framework, we construct a manifold which we call 
the lattice based multispace and which generalises, in some sense, the curve based 
multispace of Olver [50]. The multispace is a generalization of the jet bundle and 
which contains the jet bundle as a submanifold. It also contains the space of lattices 
as an open subset. The main problem with the definition of the lattice multispace is 
the fact that multivariate interpolation is not well-defined in general. To avoid this 
problem we restrict the lattices to sets of points covered by the general construction 
of deBoor and Ron ([5], [6], [7]), to what we call corner lattices. A corner lattice 
is one with just enough data to guarantee the approximation of a smooth jet. We 
restrict as well the types of coalescencing that can take place to be those along 
hyperplanes. We show that deBoor and Ron’s interpolating family is well defined 
on corner lattices and is smooth under coalescing, smooth in the sense that the 
associated Lagrange polynomials converge to the Hermite ones as the vertices of 
the lattice coalesce. Once the choice of lattice and conditions on coalescing are 
settled, we can use the interpolating coefficients to define the coordinate system in 
the multispace manifold. We notice that one can possibly consider other forms of 
lattices and coalescing, and that our theory will hold true as far as the smoothness 
of deBoor-Ron’s family is preserved. 


(i) Lagrange (ii) Hermite (iii) Taylor 
Interpolation Interpolation Approximation 


FIGURE 1. Under coalescence of the points at which the interpola- 
tion is calculated, Lagrange interpolation becomes Hermite inter- 
polation, ending with the Taylor approximation to a surface when 
all the interpolation points coalesce. By taking coordinates for the 
lattice based multispace to be the grid points and the Lagrange 
interpolation coefficients, the jet bundle is naturally embedded. 


Once the multispace is proved to be a smooth manifold, we can naturally de- 
fine a group action on it, and hence we can talk about smooth moving frames on 
multispaces. A moving frame on the lattice based multispace is, simultaneously, a 
smooth moving frame defined on the jet bundle, and a frame defined on local differ- 
ence approximations to the derivatives, depending on what point of the multispace 
the moving frame is evaluated. By defining a moving frame on multispaces, one 
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has simultaneously the full power of both the smooth and the discrete frames, and 
the smoothness of the multispace frame will ensure that we can move freely be- 
tween discrete, discrete/differential and smooth frames, ensuring that the discrete 
frame converges to the continuous one as the points in the corner lattices coa- 
lesce to create the jet. We also show that the continuity carries over to invariants 
and syzygies as well. Therefore, any smooth geometric construction carried out 
with a multispace lattice, invariants and syzygies, ensures that the final discrete, 
or discrete/differential result is an approximation of the corresponding continuous 
construction. 

In section 5 we use the multispace construction in two different applications. The 
first application is to a class of finite difference variational shallow water systems, 
which have both the correct continuum limit as well as the necessary symmetries 
for Noether’s theorem to yield conservation laws for energy, and linear and angular 
momenta, in both the finite difference case and the smooth limit. This is motivated 
by the desire to achieve an analogue of the conservation of potential vorticity in a 
numerical approximation to these equations. 

The second application concerns discretizations of completely integrable systems. 
Most well known completely integrable PDEs are linked to some geometric back- 
ground and the PDE can be interpreted as, for example, the equation induced on 
invariants by a geometric evolution of curves, or like the Codazzi-Mainardi equa- 
tions, are associated to the geometry of some type of surface. Discrete lattice 
systems also have similar interpretations ([2]). The question to ponder is whether 
or not the same geometric construction performed in the continuous case to gener- 
ate the PDEs can be carried out in the discrete case, while guaranteeing that the 
result will be a discretization of the PDE; this might be useful as a base to study 
the more interesting questions of when the discretization will be also completely 
integrable. Here we show two such processes: the Boussinesq equation is induced 
on centro-affine invariants by an evolution of star-shaped curves. We construct the 
multispace version of the construction to obtain a geometric discretization. We 
show that a modification of the construction generates an integrable discretization 
which appeared in [46]. The study of how these different discretizations might be 
related is underway. In the second example we describe the multispace version of 
the well-known construction of the Sine-Gordon equation as the Codazzi—Mainardi 
equations for Euclidean surfaces of negative constant curvature. This interpretation 
has been widely used to study pseudo-spherical surfaces as generated by solutions 
of Sine-Gordon, see [11] and [55]. Every step of the construction is guaranteed to 
discretize the continuous version, while preserving the geometric meaning of the 
elements involved. The Sine-Gordon is in fact one of several equations describing 
the surface, but which decouples from the others. In this discretization the equa- 
tions remains coupled and its integrability is not clear, but the construction itself 
is a non-trivial example of the use of mixed discrete-smooth moving frames. The 
connection between multispace and integrability is under study. 

Finally, in an appendix, we discuss a more general result concerning the dis- 
cretisation of smooth moving frames, and the continuum limit of equicontinuous 
families of discrete moving frames, with an example. 
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2. BACKGROUND 


2.1. Moving frames. Given a Lie group G acting on a manifold M with a left 
action, so that 


Gx MM, h-(g-z) = (hg): z, 


one can define a left (resp. right) group-based moving frame as a map which is 
equivariant with respect to the action on M and the left (resp. inverse right) 
action of G on itself, specifically, 


p:M—>G, p(g-z)=gp(z) (resp. p(g- z) = p(z)g™*). 


We call such an equivariant map a left (resp. a right) moving frame. The inverse 
of a left moving frame is a right one, and vice versa. 

Given a group G acting on a manifold M, the existence of a moving frame on 
the open subset U C M is guaranteed if: 


(i) the orbits of the group action all have the same dimension and foliate U, 
(ii) there is a transverse cross-section K to the orbits such that for each orbit 
O, the intersection O N K contains a single point, and 
(iii) the group element taking z € O(z) (where O(z) is the orbit through z) to 
O(z) NK, is unique. 
In this case, a right moving frame p : U > G is given by p(z)-z € K, that 
is, p(z) is the unique element of G taking z to the unique element of K N O(z). 
Since K is transverse to the orbits, the frame defines local co-ordinates given by 
ze (p(z), p(z)-z) eG Xx K. 

In the continuous case of moving frames, the manifold M could be the jet space 
J©(R?,M). In this case it is known ({17]) that provided the action is locally 
effective on subsets, as £ grows the prolonged action of G on J‘ (RP, M) becomes 
locally free. The work of Boutin ([8]) discusses what happens for products M*4 
as q grows, with G acting with the diagonal action. In any event, we make the 
assumption that for large enough dimension, there is a neighborhood of the identity 
in the group in which a moving frame can be obtained locally. 

A common way to obtain the moving frame is through a normalization process. 
One can describe normalization equations as those defining the transverse section, 
K, to the orbits of the group. If the normalisation equations are given as {® = 0}, 
then the conditions above for the existence of a moving frame are the conditions 
under which the implicit function theorem can be applied to solve ®(g- z) = 0 for 
g = p(z). Since both g = p(h- z) and g = p(z)h! solve ®(g- (h - z)) = 0, and 
the implicit function guarantees a unique solution, then p(h - z) = p(z)h~1, that 
is, p is equivariant. Typically, the normalisation equations, for which K is the zero 
set, are algebraic. Indeed, in many applications, the cross-section is a co-ordinate 
plane, so that the normalisation equations involve certain co-ordinates being set 
to a constant. Since there will be many transverse cross-sections to the orbits, 
the choice of K can greatly decrease (or increase) the calculations involved. Part 
of the “art” of the moving frame in applications is the choice of cross-section, or 
equivalently, the choice of normalisation equations. 

Given a moving frame (left or right) one can generate all possible invariants of 
the action. Indeed, if p is a right moving frame, the expressions 


plu) -v 
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for any u,v € M are clearly invariant; their coordinates are called the basic 
invariants. One can easily see that any invariant of the action is a function of 
these, using the replacement rule: If I : M — R is invariant under the action, so 
that I(g-v) = I(v) for all g € G, then setting g = p(u), one obtains 


I(p(u) + v) = Iv). 

Different choices of the manifold M gives rise to different familiar cases. For ex- 
ample, if M is the jet space J‘) (R?, P) for some manifold P where G acts, and G 
acts on M via the natural prolonged action given by the chain rule, then p would 
generate moving frames on p-submanifolds and the invariants will be standard dif- 
ferential invariants (for example, curvatures, torsions, etc). If M = x,P is the 
cartesian product of a manifold P where G acts, and G acts on M through the 
diagonal action, then the invariants are the so-called joint invariants (see [51]). 

The authors of [41] defined discrete moving frames, essentially a choice of group 
element associated to each vertex in an equivariant way. The discrete moving 
frame can be defined to act naturally under the shift operator, greatly simplifying 
calculations with discrete (difference) equations. We next review this definition. 

Let GN denote the Cartesian product of N copies of the group G. Allow G to 
act on the left on G using the diagonal action g - (gr) = (ggr). We also consider 
what we have called the “right inverse action” g - (gr) = (grg~'). 


Definition 2.1 (Discrete moving frame). We say a map 

p: MN > GN 
is a left (resp. right) discrete moving frame if p is equivariant with respect to the 
diagonal action of G on M% and the left (resp. right inverse) diagonal action of 


G on GN. Since p((a,)) € GY, we will denote by ps its sth component, that is 
p = (ps), where ps((£r)) € G for all s. Equivariance means, 


P(g > (2r)) = as((9- tr)) = gos((tr)) (resp. p5((xr))g7") 


for every s. Clearly, if p = (ps) is a left moving frame, then P = (p7!) is a right 
moving frame. 


As in the original group-based moving frame definition, if (us) € M, one can 

define 
Ps Up =I. 

The coordinates of these invariants for any r generate all other invariants even 
when s is fixed (see [41]). We note that the action induces an action on the co- 
ordinate functions, the same as it induces an action on any function, specifically, 
g: f(u,) = f(g-u,). The components of IT will be invariant as IT is, they are called 
the basic invariants. 

One can describe a smaller set of generating invariants, the Maurer—Cartan in- 
variants. 


Definition 2.2. Let (ps) be a left (resp. right) discrete moving frame evaluated 
along a twisted N-gon, the element of the group 
Ks = p;'ps+1 (TeSp. penis”) 


is called the left (resp. right) s-Maurer-Cartan matrix for p. We call the equation 
Ps+1 = Ps Ks the discrete left s-Frenet-Serret equation. 
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If G C GL(n, R), then the entries of the Maurer—Cartan matrices, together with 
the basic invariants I$, generate all other invariants. See [41] for more details. 


2.2. Lattices. Lattices are subsets of Z? with a variety of properties. We first 
define adjacency. 


Definition 2.3. Two points m, n € ZP are said to be adjacent if 
J 


Definition 2.4. We say that a subset [ C Z? is a connected lattice if it consists of 
a single point, or, if between any two points mı, mz € I there is a path, m = nj, 


Nj, ... Ny = Mz, such that n; is adjacent to nji; for i =0,...,N—-1. 
(i) 4 (ii) + 
a- + + a- 
+ + a- a- H- 
a- Ar 
cii) (vw) E 
+ a a 
+ 
at- + a + + 
_ a 
+ 
++} +—. 


FIGURE 2. In (ii) and (iii), T is connected. The subsets (i), (iv) 
and (v) are not connected lattices as there are points not adjacent 
to any other point. 


Natural operators on Z? are the well-known shift operators, namely {7;}?_,, 
where 
Tig sks NG i i = (Nizi gig Pica g Mg 


We will also consider lattices in a manifold. 


Definition 2.5. The image of a lattice IT C Z” in a manifold M by a map ®:T > 
M is denoted by L, and is also called a lattice. We assume this map to be injective, 
a condition which will be relaxed under controlled conditions in the second part of 
this paper. The adjacency in L is determined from that of T. 


We will assume from the beginning that the lattice L C M is contained within 
one coordinate chart for the manifold, so that for all practical purposes we can 
assume the manifold is R”, or often a parametrized surface with the parameter 
serving as local coordinates. We will also assume that lattices are connected (one 
can study each connected components separately). 
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The lattice itself does not need to be covered with one lattice neighborhood, 
however. We define next a lattice variety, which will allow us to work on lattice 
models of spheres and tori. 


Definition 2.6. A lattice variety £ is a set of points that can be covered by subsets 
La, each of which mapped by an injection a to a connected lattice La C ZP 
for some fixed p. Every adjacency in £ is contained in at least one of the La. 
Furthermore, in the overlap La N £g, the gluing map ġa ° os preserves adjacency. 
We call (La, Ga) a local lattice coordinate system. 


Lattice coordinates essentially introduces an order in the lattice (inherited from 
ZP through ġa) so one can clearly define shifts. We say that a shift map is defined 
at a point in £, if it is defined in at least one chart. Since the chart interchange 
maps preserve adjacency, the existence of a shift map is well-defined. 


3. MOVING FRAMES ON LATTICES AND LATTICE VARIETIES 


Let £ be a lattice variety and let N = |I| be the number of vertices in £. Let 
Ly be the set of p-lattice varieties in M with N vertices (in much of what follows 
N could be infinite). 


3.1. Moving frames, invariants and Maurer—Cartan invariants. Let G bea 
group acting on Ly (for example, if the lattice lives inside a manifold with a group 
action, the action would be the one induced on the lattice), and for simplicity 
assume that it is a left action (that is, g-(h- u) = (gh)-u. A parallel description 
can be made for right actions). 

A discrete moving frame will associate an element of the group to each vertex in 
the lattice in an equivariant fashion. 


Definition 3.1. [Moving frames on lattices] 
Let U be a subset of Ly. We say 


p:u—+GN 


is a left (resp. right) discrete moving frame on U whenever p is equivariant with 
respect to the action of G on £y and the left (resp. inverse right) diagonal action 
of G on GY. That is, if p = (p;)§_, denote the components of p in G^, then 


pilg-L) =gpi(L) (resp. plg- pi(L)) = pi(L)g~*), 
with |T| = N. 


In general, moving frames exist only locally, which is the reason why we need 
to restrict its domain in Ly. Given a lattice variety £, and a coordinate system 
indexed by I C Z’, the moving frame p assigns group elements at each vertex (thus 
N of them). We will call pr the moving frame at the vertex R €T. Note when the 
index is applied we are assuming the use of local variety coordinates. 

It is a simple matter to go from a moving frame to a discrete moving frame by 
taking a family of cross-sections, one per lattice vertex, as stated in the following 
result. 


Proposition 3.2. Let {Sr| RET}, be a family of sections, with Sp transverse to 
the orbit of G at up, assumed to be indexed locally by T. Let g = (gr) € GN be 
uniquely determined by the condition 


(1) JR'LESR 


DISCRETE MOVING FRAMES 9 


for L € Ly and R ET in some coordinate system for L. Then (gr) = (pr) is a 
local right discrete moving frame. 


The proof of this statement is straightforward from the discussion of the moving 
frame. We note that a moving frame is defined as an element of the group on the 
entire lattice, but normalization equations give frames defined only locally. The 
interchange maps from one domain to another will play a role in what follows. 

Remark on Notation. From now on a multi-index will denote the use of local 
lattice coordinates, while the lack of it will indicate global definitions. Also, we will 
denote by pr each individual component of p, or the moving frame at the vertex 
UR. 

Example 3.3. Consider 1-lattices - or polygons - in the Euclidean plane. The 
group E(2) can be identified with the subgroup of GL(3, R) given by 


2) jŠ G a) 


with © € O(2) and b € R?. It acts on R? as 


(0) (2)= (ers) 


with x € R?. We choose as our transverse cross-section, the one given in coordinates 
by Pn-Un = 0 and py: Un41 = €1|Aun||, where Aun = Undi — Un. If pn is displayed 
in the matrix representation as in (2), solving the system 


Pn ` Un = Onun T bn = 0, Pn ` Un+i = OnUn+1 + bn i 0, Aun = Aun jer 


results in the left moving frame 


An 1 0 ) o ( 1 0 0 ) 
Pn \-O71bn Oz! |Aunler aat -IAT 


where J = & is the canonical symplectic matrix. 


Definition 3.4 (Invariants and basic invariants). We say the function on Ly 
Ds Ln >R 
is a lattice invariant under the action of G if I(g-£L) = I(L) for any g € G, LE Ly. 
A local invariant will have the same property in some coordinate chart. 
Given a right moving frame p on p-lattices, we call the invariants 
Ir =pr(L)-L 
the basic invariants, where p(L) = (pr(L)) with pr(L) € G (the group element at 
the vertex in position R depends, in general, on all vertices of the lattice). Once 
we choose coordinates in £, given by (uz), the local invariants are defined to be 


IR = pr((us))- um 
for R, M eT. These are clearly invariants of the action. 


The basic invariants take their names from the fact that they generate all other 
local invariants. In fact, they do generate them even when R 1s fixed. 


Proposition 3.5. If I is any lattice invariant, then I can be written as a function 
of the basic invariants IN for any fired RET. 
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Proof. This is an immediate consequence of so-called replacement rule. Let (un) 
represent the vertices of the lattice variety. If I = I((uy)) is an invariant of the 
action, then I(g - (un)) = I((un)) for all g € G, and in particular for pr((u,z)). 
Thus 

I(pr((us)) > (un)) = IOR) = L((un)) 


which shows us how to write J in terms of the basic invariants with R fixed. 


From this wealth of invariants we will be selecting a few, the so-called Maurer- 
Cartan invariants. They are the discrete analogue of the invariants defining the 
classical Frenet—Serret equations and, like their continuous counterpart, they will 
form a generating system. 

From now on we will extend the shift operators in the standard way to algebraic 
functions of wy using the properties 7;(uyur) = TiusTiur. We can also apply a 
shift to pz by sending pJ to py+e, or we could apply it to py by shifting the variables 
ur that pz depends on. But notice that unless the sections in 1 are shifts of each 
other (i.e. Spye, = TiSr for all i = 1,...p and all R € T), these two operations 
do not need to produce the same result. Given that in many of our situations and 
in all of our examples we do assume the sections to be invariant under the shifts, 
we will abuse the notation and denote all these maps by 7;, so that, for example, 


TiPs = Pi+e;: 


Definition 3.6. [Maurer—Cartan invariants] Let p be a right moving frame along 
p-lattices. We define the right (R,i)-Maurer—Cartan group element to be KR,i), 
the element of the group given by 

Kiri) = (TIPR) PR = PR+e:PR - 


Its local coordinates (or the entries of the matrix, if G C GL(n,R)), will be called 
the (R,i)-Maurer—Cartan invariants. 


Theorem 3.7. Let p be any right moving frame. The (R,i) Maurer—Cartan in- 
variants, i=1,...,p, RET, together with pr((us))-ur = IB, RET, generate all 
other invariants for the action of G on Lyn. 


Proof. The proof is based on what are commonly known as the recursion formulae. 
Directly from the definitions we get that 


Kray Ih = (Tier) bp’: (Pr: Um) = Pree; UM = IRh e, 
and from 
(3) Kaa Ip = Ih 


we have 


M —1 M 
If = Kay Th 


Now, since T is connected, given M € ZP, any R € ZP is related to M € ZP through 
either recurrently increasing or decreasing its individual components, using the shift 
operator. At each step the invariant obtained when increasing or decreasing the 
components in M is generated by those in previous steps and by Maurer-Cartan 
invariants. Thus, we can start using IM and reach I M , for any R, using both 
versions of the recursion formulas. This proves the statement of the theorem. 


te; 
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Example 3.8. The simplest example is the translation group viewed as a subgroup 
of GL(n+1,R) and acting on R” as 


peall D UAT er) mete 


where n is the projection in the last n components. If (ur) is a p-lattice in R”, 
a transverse section to the orbit of the group at ur is given by ur = 0. Thus, 
the moving frame is determined by g: ur = 0, which implies a = —ur. The 
basic invariants are ie = PR: Uy = Uy — Up, while the Maurer—Cartan matrices 
are PRH PR whose only non-constant entries are the Maurer—Cartan invariants 
UR— Ure, RET, i=1,...,n. Itis straightforward to show that the IM can be 
written in terms of the Maurer-Cartan invariants Les, 

Example 3.3 cont. In the case of the Euclidean plane, we found a moving right 
moving frame given by 


(4) TE E D 


where 


—1 Un Un 
O, = (aa Ja) 


The basic invariants in this case are given by 


for any n,m. Notice that pp -Un41 = |Aunle:. The Maurer-Cartan matriz is given 
by 


(5) 1 1 0 
Pn+1Pn = —On419bn + bn41 On+197! 
where 
1 —u - Au 
-On Onbn bn = -On a — un) = n+1 n 
+1 + On+1 +1(Un41 — Un) [Aural Ges) 
and 
O ETEA N 1 Un+1° Un det (un, Un41) _ ( cosan sin Qn 
mttin © [Atnsil Aun] \- det(unti un) Unga tn -sinan cosan)” 


Where ay, is the angle between un+ı and un. Therefore, a generating set for the 
Maurer-Cartan invariants are |Aun|| and an, for alln. Since the basic invariants 
are also generated by |Au,|| and the angle between un and um, the Maurer-Cartan 
invariants generate all the basic ones. 


3.2. Syzygies. 


3.2.1. Basic local syzygies. From the definition of Maurer—Cartan matrix, KN) = 
(Tien) py’, we have 

TK 3.) = GRD = (TTips) er Koz) 
and also 

TK (3,3) = (Tips) Tips = (Tiener Kay 
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Given that shifts commute, we obtain 

(6) (HK) Kay = (TKa) Kua. 

This expression gives us a number of algebraic relationships between the different 
Maurer—Cartan invariants. We will refer to these as basic local syzygies, a discrete 
generalization of the differential syzygies that are satisfied by differential invariants, 
such as the Codazzi—Mainardi equations for Euclidean invariants defined on sur- 
faces. The above syzygies generate most of the possible algebraic relations among 
Maurer—Cartan invariants. Further independent syzygies may be created by the 
topology of the lattice L. 


Definition 3.9. We say two syzygies are equivalent up to basic syzygies if one of 
them is an consequence of the other and syzygies of the form (6). 


To describe global syzygies, we need first to define the discrete fundamental 
group of the lattice and related material standard in the study of the topology of 
graphs [47]. 


3.2.2. Discrete fundamental group. Let L be a p-lattice in M, that is, the image 
of a map from a connected p-dimensional subset I of ZP to M. As before, we will 
denote the image of J € T as (J) = uy € M. For simplicity, let us assume that 
T = Z’, although one can apply much of what we will say next to other cases. 
Notice that, in principle, we are allowing cases when the map ¢ is not 1-to-1 so L 
does not need to have a trivial topology. 


Definition 3.10 (Paths). We say a subset P C L is a path of length r joining two 
points a,b € M, if it can be ordered as P = {z;}'_9, with £o = a, £r = b and 2; 
adjacent to x;41, for all i = 0,...r — 1. We say the path is closed if a = b; we say 
it is simple if x; A £; for any i # j (except perhaps £o = 2, if closed). Notice that 
by giving the vertices of the path in a certain order we are implicitly assigning an 
orientation to it. This will be relevant once we associate syzygies to closed paths. 


Definition 3.11 (Edge). Given a path y in a lattice, y = {x;}7_), we say the 
ordered pair [x;,2;41] is an edge of the path. The ordering gives an orientation of 
the edge. 


Definition 3.12 (Sum of paths). Consider the set of all closed paths with base 
point a. One can define the sum of two such paths by concatenation; that is, if 
{zi}f-9 and {y;};_9 are two paths, their sum is given by {xi} + {yj} = Sopher 
with 
Zi = Ti, Tes Zr+j+1 = Yj, J= 0 yc ony 
If the paths are not closed, but x, = yo, one can equally define the sum of the paths 
by concatenation. 
Definition 3.13 (Basic homotopy). A transformation of a path {x;}, ¢({a:}) = 
{y;} is a basic homotopy if {y;} is equal to {x;} except for 
(1) adding or removing a subpath of the form [2,244] + [@r41, £r]; 
(2) changing a subpath of the form [z,., £r+1] + [£r+1, £r+2] by one of the form 
(tr, Yr+3] + [yr+3, £r+2], where 2,41, Ur+2, Yr+3 form a basic square of 
the lattice. 


Transformations (1)-(2) are called the two basic homotopies. 


Figure 3 shows examples of basic homotopies. 
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FIGURE 3. Examples of the two basic homotopies applied to the 


blue paths 
The lattice is missing the inner vertices: the red path 
Decomposition of two homotopies into wraps around the missing tube while the blue one 
basic homotopies is homotically equivalent to the constant path at a. 


FIGURE 4. The red and blue paths on the left are homotopically 
equivalent, with the yellow and green paths showing the sequence 
of basic homotopies required, while the ones on the right are not. 


Definition 3.14. We say the two paths {x;}/_) and {yj}; joining a and b, are 
homotopically equivalent if £o = a = yo, £r = b = ys and either the paths are equal, 
or one can be transformed to the other by a finite sequence of basic homotopies. 


By construction, homotopy of paths joining a and b is an equivalence relation. 
Figure 4 shows pairs of homotopic and nonhomotopic paths. 


Definition 3.15 (Discrete fundamental group). Consider the space of all closed 
paths based at a, and let 7,(LZ,a) be the set of homotopy classes of these paths. 
The operation above endows 71(L,a) with a group structure. We call mı(L, a) the 
discrete fundamental group of L. 


The fundamental group does not depend on the point a chosen, as far as the 
lattice is connected. (It suffices to join a to a different point b using a path P, and 
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use P to relate closed paths based on a to those based on b by conjugation, as 
done in the continuous case.) 


3.2.3. Syzygies associated to closed paths on a lattice and global syzygies. Assume 
we have a moving frame along a path. To each edge of the path we can associate 
a Maurer-Cartan matrix of invariants as follows: 

Assume either that x; = uy and 2441 = Thuy = Us+e,, or that £ti+ı = T; ‘uyg = 
us—e,- To [%;, 2141] we associate the matrix 


(7) Korn) = Te(ps)P7) = Pr+enP 7 


in the first case, and 


= _ -1 _ 
(8) Kj enk) = (Tel PJ-er)P7L en) = beeps 


in the second case. 

We note that the choice of Maurer—Cartan matrix depends on the orientation of 
the path, with an edge being associated to the inverse matrix if the orientation is 
reversed. 

Next we will associate a group element K(y) to each path y on the lattice, 
namely, the product of the Maurer-Cartan matrices along the path. From the 
definition of the Maurer-Cartan matrices in terms of the discrete frame, it will be 
evident this product telescopes to involve only the discrete frame at the endpoints of 
the path, provided the discrete frame is defined along the whole path. Nevertheless, 
in terms of the components of the Maurer-Cartan matrices, the Maurer-Cartan 
invariants, the product will not telescope — this is the syzygy. Evaluating K (y) for 
closed paths y, leads to relations on the invariants. It becomes important to find 
those relations which are non-trivial, in the sense that they are not an algebraic 
consequence of basic local syzygies, given in Equation 6. 

In what follows, we will show that K(y) is a homotopy invariant. Evaluating 
K(y) on closed paths which are not homotopic to the constant trivial path leads 
to relations on the Mauer-Cartan invariants which cannot be obtained in terms of 
the basic local syzygies. 

We start with paths which lie in the domain of a discrete frame. This, of course, 
need not be the case since the existence of moving frames is guaranteed only locally. 
If we need to cover the lattice with several co-ordinate patches on which discrete 
moving frames exist, we will obtain invariant transition matrices associated with 
the cover of the lattice defined by the domains of the discrete moving frames. We 
discuss this more involved case later in this section. 


Definition 3.16. Assume a discrete frame exists along a path y. The product 
of the Maurer-Cartan matrices along y is denoted K(y). Specifically, for the path 


{xi}, we have 


K ({xi}p2o) = K([ar-1, £r]) + K([20,)), 
where K((x;, £i+1]) is the Maurer-Cartan element associated to the edge [;, £i+1] 
as in (7)-(8) so that K([a;,vi41]) = K7*[vi41,2;]. If y is the constant (trivial) 
path, we define K(y) =e. 


It is evident that if y = y, + %2, then K(y) = K(y2)-K(m1). (See Definition 
3.12 ). 
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= a 
K2=P)+eķ+ejPJ+ek 


X2=Ujie, X2=U 
$r — te ^3 M J+ej+ek 
A 
K3= a 
1 3 PJ+ej PJ+ej+eķ 
Ky=Pj+e,Py 
X,=Uy X4=U+ej 


K4=P) Piso 


FIGURE 5. In this Figure, we see that along the closed path 
(@1,%2,3,04 = £1) we have K,K3K2K, = e, the identity in G. 
This is equivalent to the basic local syzygy, Eqn. (6). 


In Figure 5, we illustrate the basic local syzygy, in the form K({x1, £2, £3, £1}) = 
e along a closed path of length four. 


Proposition 3.17. Let yı and y2 be two paths joining a and b in the lattice L 
for which a discrete moving frame exists. Assume that yı is homotopic to y2, then 
K(m) = K(2). 


Proof. Since a homotopy is a finite composition of basic homotopies, it suffices to 
show that if y, and y2 differ by a basic homotopy, then K(71) = K (72). 

In the case (1), this is trivial since the only difference between K (q1) and K(72) 
is a product of the form K;K;' = e. In the case (2) it is equally simple since 
they differ only by a product K;Kj41 appearing in K(y,) and K;Kj+1 appearing 
in K( 72), with K;Ki+ı = K;Kj41 being a local syzygy since their vertices form a 
square in the lattice. 


The following Corollary is an immediate consequence of the previous proposition. 


Corollary 3.18. If our lattice variety is covered by one coordinate system and 
there exists a moving frame defined everywhere, each class |y] of the fundamental 
group of the lattice defines what we will call a global syzygy of the Maurer-Cartan 
invariants, in the form K(y) =e. 


For example, if a discrete frame is defined on an annular lattice L, then there 
will be a path not homotopic to the constant path within L, but K(y) = e for all 
closed paths. 

We now turn our attention to paths that move through different domains of 
discrete frames. We can still associate a group element K to the path by patching 
local products of Maurer-Cartan matrices, assuming the local neighborhoods where 
the different discrete frames are defined overlap. Assume we have two subpaths 
{@o,.-., vi} and {x;,vi41,...,2;} of a closed path, and assume we can find a 
moving frame ps at zs, s = 0,...įi and a different moving frame fs for s = i,i + 
1,...,j (notice the overlap in the subindices). We define M(x;) = f;p;' so that 
ĝi = M(a;)p;. Clearly, since both p; and ĝ; are right equivariant, the matrix M (x;) 
is invariant. Then, to the path y = {xo,...,x;} we can associate the product of 
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invariant matrices 
K(y) = (6j0j"1)(Aj-10j 2) --- (Bit; ")M (z2:) (Pip) - -- (0100) 
= Ky-1...KiM(a)Ki-1 ... Ko 


with the invariant matrix M (x;) linking the Maurer-Cartan matrix in one coordi- 
nate system to the next. 


Definition 3.19. If a vertex z lies in the domains of both the discrete frame p and 
the discrete frame p, we say the group element 


(9) M(x; P, p) = p(x)p(x)~* 
is the transition Maurer-Cartan matrix at the vertex x, associated with the change 
of frame from p to P. 


Important Assumptions. We already have the assumption on our co-ordinate 
charts La, that every edge appears in at least one La, and so every Maurer-Cartan 
matrix can be written in (at least one) co-ordinate system. We assume further that 
every edge is in a domain of a discrete frame. In this way, every Maurer-Cartan 
matrix, every transition matrix, and their local products, can be expressed with 
respect to a single set of co-ordinates. By taking a refinement of our co-ordinate 
cover as necessary, we therefore assume that our cover consists of sets which are 
domains of both frames and co-ordinate charts, and that every edge appears in at 
least one element of the cover. 

If a vertex lies in the intersection of several frame domains, there will be co-cycle 
conditions. For example, if x € dom(pa) N dom(pg) N dom(p;5), then clearly 


(10) M(x; ps, pp)M (2; pp, Pa)M (a; pa, ps) = e. 


On the other hand, two equal paths can have different group element K(y) if 
the choices of either the initial or the final moving frames are different. Thus, our 
element of the group does not only depend on y but also on the initial and final 
choice of moving frame. In this case we will denote the group element above 


K(y, a, 8) 


for the element of the group that starts in dom(pq) and ends in dom(pg), and 
where these are the choices of frame for the calculation of the initial and final 
Maurer-Cartan matrices. By analogy, we will also denote by M(a;a,) the 
matrix M (z; pa, pg). 


Lemma 3.20. K(y,a, 8) does not depend on the choice of coordinates or moving 
frames one chooses along y, only on the initial and final ones. 


Proof. Assume at some point x we make different choices of moving frame at a 
point x in the overlap of (at least two) different frame domains, so we move from 
Pu tO py for one path and from p, to p, for the other. In that case we introduce the 
transition factor M (x; n, 4) in one of the lifts, and M (x; v, u) in the other one, and 
we continue the different paths using the corresponding choices. At some point we 
need to come back to a common choice, even if that happens only at the end of the 
path. But when we expand the different factors of K(y,a, 8) in terms of moving 
frames, the intermediate factors all vanish as we saw before, and the difference 
is only at the beginning and at the end of the product. Thus, without losing 
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\ / Yo 


š PAN 
se” Domain for f 


FIGURE 6. Several ways a path can move through an interchange. 


generality we can assume that we come back to a common moving frame right after 
we introduce the split. 
That is, a path includes the factor 


M(x; 8,)M(x;n, 1), 
while the other includes 

M(x; 8,v)M (a; v, u). 
But using (10) we have that both these factors are equal to M(x; 8, u), and hence 
they are equal. 


We now argue that K(y,a, 8) is still a homotopy invariant for paths that start 
and end with the a and 8 choice of moving frame, even if the path moves through 
changing domains of discrete frames (and changing co-ordinate systems). Consider 
Figure 6, in which we assume that [x,y], [y, z] are edges in the domain of p, while 
[x,t] and [t, z] are edges in the domain of p. We consider paths q; running from zo, 
in the domain of p, to yo in the domain of P. If yı goes via x and t, then we must 
switch from p to p at x. We achieve this by using the matrix M(x) = p(x)p(x)~! 
(which is invariant as the frames are right frames), and then defining K (71) to be 


(lt, yo])) E ([x, #]) M (x) K ([xo, 2]) 
yo)alt)* (PAPE) (A(x) p(x)~*) p(a)p(«o)~! 
yo)p(xo)*. 


Considering the path yz from zo to yo via z results in 


Ky) 


R N) 


>) 


nN nN 


K (72) = K (|t, yol) K ([z, t]) M (2)K (Ly, 21) K (xo, y]) 
which is also equal to A(yo)e(xo)~* by a similar argument, and thus we have 


(11) K (1) = K2) 
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FIGURE 7. The situation of Figure 6 translated to £, the “lift” of £. 


Since we have homotopy invariance of the K element within domains, in this way 
we can see that even passing through a change of domain, we maintain homotopy 
invariance. 

In order to prove homotopy invariance in general, we construct a “lift” L, of 
L, with respect to a good cover of L, using which we can keep track of which 
discrete frame we are using on our paths in £. The lift £ that we construct is 
not a lattice variety in general, and does not lie in M, but nevertheless serves our 
purpose here. To construct L, we take the disjoint union of the charts, together 
with their edges (adjacencies), and for every  € La N Lg, we take a new adjacency, 
or edge, [æ|a, x|], and let this be a new edge in Z, with the associated Maurer- 
Cartan element being pgp,', the transition Maurer-Cartan matrix. See Figure 7, 
which shows the lifting for the case of Figure 6. Define the projection from £ to £ 
as the natural projection that collapses the different copies of the vertices which 
lie in intersections of charts. That is 


ihe. aoe E if x; belongs only to one domain 


x if a; =x, where zo is a lift of rE La 


A path 7 in £ is a lift of a path in £ when the projection of y is y. Lifts may 
not be unique, as they depend, for example, on the element of 7~!(a) at which 
the path begins, where a is the initial vertex of y. 

Changing from one frame to another along a path is, in Z, simply proceeding 
from one vertex to another, with the transition Maurer—Cartan elements, Definition 
3.19, being the group element associated to the new edge. We note that a change of 
co-ordinates simply changes the local labelling of the points, and so is less important 
when considering the group element associated to a path. 

We now show that in £ we have new basic local syzygies obtained from the 
transition Maurer-Cartan matrices. These new local syzygies allow us to obtain a 
result similar to that of Proposition 3.17, but for the more general situation where 
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path t,x,y,z in L Lift of path in L 


FIGURE 8. New edges in £ give rise to new local syzygies, here, 
K (lye, £6])K (|e, £a]) = K (lyg, Ya])K (lya, tal) 


we need to change domains of our moving frames. Consider Figure 8. Let the 
domain of p be La and the domain of p be £g. Consider the simple closed path 
1 = {zla £p, Yle: Yla Tla}. If K (lys, x6]) = ply)e(2)~*, K (las, tal) = Alx)p(x)™, 
K([ys,Yal) = Py)p(y)~*; K (lya, tal) = ply)e(z)~*, the K element for this path 


is, 
K(7) = (PE) P) PPA) PAE) (Az) p(a)~") = e 
showing this path defines a basic syzygy in £. 

Finally, we define the monodromy of a closed path in £. When such a path is 
lifted to £ it need not be closed, as it may begin in one frame domain and return 
in another. Consider the lift 7 of a path y beginning at xo in the domain of the 
frame pa, and ending at xo in the domain of the frame pg. We can close the path 
in £, by adding to 7, the edge [zg,rq], but in principle this need not happen. 
Therefore, the element of the group associated to the lift would not be equal to e, 
but rather to M (x; a, 8) (or its inverse, depending on the orientation). We call this 
group element the monodromy of the lift. 

We can now state the more general Theorem concerning homotopy invariance. 


Theorem 3.21. To each path y in the lattice variety L, and to each choice of 
initial and final moving frame, we can associate a group element, K(y,a,(), such 
that if yı and y2 are homotopic in L, then K(%1, a, 8) = K(y2,a, B). 


Proof. We first note that £ has essentially a global moving frame, in the sense that 
every vertex has an equivariant group element associated, namely pa for £a and pg 
for xg, and so on. Hence we can define K (y) for any path, the same way we did 
previously. 

We first note that any lift of the trivial path to a closed path in £ has, by the 
co-cycle conditions, Equation (10), K(@) = e. When we lift a path y to £, we need 
to keep track of which discrete frame we are working in, but as we saw before, the 
element K(7,a, 8) is only affected by the beginning and end choices. 

Assume two paths are homotopically equivalent, and let’s lift the homotopy. 
By construction of £, we can assume the endpoints of the lift of the homotopic 
paths also remain fixed and determined by the a and 8 choices. As in the proof of 
Proposition 3.17, it will suffice if we show that two paths that differ by one of the 
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basic homotopies have the same K(7,a, 3), even if we need to change the moving 
frame domain. But this was already proved in the argument concerning Equation 
(11). 


The syzygy of any closed path y, where the domain of the discrete frames are 
considered to be the same at the end and at the beginning, is K(7,a,a) = K(y) = 
e, while those where a # 6 will have a non trivial monodromy K(y,a,() = M #e. 
Furthermore, different choices of a@ give rise to the same syzygy: the group elements 
are related by conjugation 


K(7,a, a) = pops K (7, 8,8) ppp a's 


and the transition matrices Pog are essentially a change of coordinates. Indeed, 
a frame defines a local coordinate system of the form U x K where U is a neigh- 
bourhood of the identity in G and K is the cross-section which has invariants for 
coordinates ([37], Chapter 4.) We thus have the following Corollary. 


Corollary 3.22. Let |y] be the homotopy class of a closed path y in L. Each 
element |y] of 7(L) gives rise to a syzygy on the Maurer-Cartan invariants, in the 
form K(y) =e. 


Notice that associating a syzygy to a closed path with a monodromy is essentially 
the same as associating a syzygy to the closed path for which the beginning and 
end moving frames are the same. Indeed, K(y;a,a@) = K(y,a,8)M(a;6,a), and 
hence K(y;a,a@) = e is the same syzygy as K (y; a, 8) = M(x; a, 8). 


Example 3.23. Consider a bi-periodic lattice L = {znm}nmez With Zn+4kTi,m = 
Zn,m and Zn, m+kTə = 2n,m for some periods Tı, Tə € Z. For simplicity, assume 
we can find a global moving frame p = (Pn,m). The topology of this lattice is 
comparable to that of a torus, and one can easily show that 


m™(L) = LZ. 


The two generators of mı(L) (marked in Figure 9 with different colors) are repre- 
sented by the two global syzygies 


Ko,mKı,m see KT, —1,m =G Kn oKn,1 tee Kn Ta—1 =œ. 


If a global moving frame does not exist, then the product might be equal to a 
monodromy matrix that will depend on the choice of moving frame at the beginning 
and the end of the closed path. If we choose the same moving frame, then the 
syzygy will be independent of the point chosen as beginning and end, and it will 
be as above. 


4. CONTINUOUS LIMITS OF THE DISCRETE PICTURE: LATTICE BASED 
MULTISPACE 


In this section we show how one can construct a continuous moving frame as 
the limit of a discrete one, and vice-versa, by coordinating the transverse sections 
that determine them. This is achieved using our lattice based multispace, in which 
derivatives and their finite difference approximations exist in a single manifold con- 
taining both the jet bundle and cartesian products of the base space. Both smooth 
and discrete frames are then part of a single frame on this multispace, and their 
relationship is given by the continuity of the multispace frame under coalescence. 
We show in this case that not only moving frames but also discrete invariants and 
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FIGURE 9. Discrete toroidal lattice variety, with periodicities T} = 
8 and To =6 


local discrete syzygies converge to differential invariants and differential syzygies 
respectively. We will use our multispace constructions to describe discretizations 
of integrable systems, and to finite difference models of variational systems, in §5 
In the Appendix, we will show that more generally, an equicontinuous family of 
discrete frames will converge to a smooth frame. 

First of all, we recall the definition of multispace in one variable, as developed 
by P.J. Olver in [50]. Olver provides co-ordinates for his space of pointed curves 
in terms of the Lagrange approximation of the curve via interpolation at specific 
(given) points. The coordination of the discrete and the smooth pictures is a conse- 
quence of the fact that Lagrange interpolation becomes Hermite interpolation under 
coalescence. In order to provide coordinates for our higher dimensional generali- 
sation of Olver’s construction, we need to restrict our ‘pointed surfaces’ to those 
where the Lagrange and Hermite interpolations are similarly related and to where 
the interpolations vary in a smooth manner with respect to the data. The details of 
the interpolation are critical, since the coefficients of the interpolation polynomial 
will define the desired coordinates. We use the theory of multivariate polynomial 
approximation due to de Boor and Ron, [5], [6] and [7], described in §4.2. We then 
describe our lattice based multispace and prove that it is a manifold. In fact, we 
detail two related versions of multispace, one containing the jet bundle J(M,R) and 
one containing the jet bundle J(U, M) where U C R? for any p < dim M. Both 
arise in the applications. Thereafter we show that the limit of the discrete Maurer 
Cartan matrices are the smooth Maurer Cartan matrices and that the local syzygies 
(Equation 6) limit to the so-called zero curvature condition of the smooth Maurer 
Cartan matrices, in §4.4. We also describe mixed discrete/continuous cases. 


4.1. Multispace in one variable. The idea behind the definition of multispace is 
to create a manifold where both discrete and continuous cases coexist in one over- 
arching smooth construction, where the continuous frame is a limit of the discrete, 
and the limit of the discrete data is the continuous data. Multispace resembles the 
jet spaces, but includes also discrete versions of the jet spaces. 

Given a manifold M, we let C™ = C™ (M) denote the set of all (n+ 1)-pointed 
curves contained in M; that is, the set of (zo, ..., Zn; C), where C is a curve and 2; 
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are n+ 1 points in C, not necessarily distinct. We denote by 
ni = #19 | zj = zi} 
the number of points that coincide with z;. 
Definition 4.1. (Multispace for curves) Two (n + 1)-pointed curves 
C=: (205 0225 2nC); C= (Zo <i Bar O), 
have nt? order multi-contact if, and only if 


= pn ila foreach i=0,...,n, 


Zi = Zi, and jn,;-1C 


where jC denotes the kth jet of the curve C, and where >, n; = n. The n* order 
multi-space, denoted M”), is the set of equivalence classes of (n+1)-pointed curves 
in M under the equivalence relation of nt order multi-contact. The equivalence 
class of an (n + 1)-pointed curve C is called its nt? order multi-jet, and is denoted 
by jnC E€ M”. 


When all the points are distinct, then two curves belong to the same equivalence 
class whenever they have the distinguished points z; in common. Thus, we can 
identify this special subset with the off-diagonal Cartesian product, denoted by 
Me"+) in [50]. On the other hand, if all the points coincide, then the class is 
equal to the jet class. Thus, both extremes can be found in one space, together 
with all the intermediate cases. In the first part of [50] the main result is the 
following theorem. 


Theorem 4.2. If M is a smooth m-dimensional manifold, then its n™ order mul- 
tispace M™ is a smooth manifold of dimension (n +1)m, which contains the off- 
diagonal part M°\"*)) of the Cartesian product space as an open, dense submani- 
fold, and the n™® order jet space J” as a smooth submanifold. 


The topology is inherited from that of the manifold M, and the proof is based 
on finding coordinate systems in a neighborhood of an equivalence class C. The 
coordinate system is given by the classical divided differences and their limits. That 
is, given a curve C with n + 1 distinguished points {z9,...,2,} and with a certain 
order of contact n; at each point, there exists a unique polynomial p of degree 
n such that p(z;) = C(z;) and such that p™ (z) = C™(z;) for any k < ni and 
any i = 0,...,n. The polynomial is a natural representative of the class C and 
its coefficients provide smooth coordinates in a neighborhood of C. Of particular 
importance is that the coordinates are smooth under the coalescence of points zi. 
For more details, see [50]. 

In the second part of the paper [50], the author assumes there is a Lie group G 
acting on the manifold M, and he defines the action of this group on the multispace 
as that naturally induced by it: the action on the differential part is the prolonged 
one, and explicit formulae for the action of the group on classical divided differences 
are given. He also explains how, assuming that one chooses a cross section to 
the orbit of the group at a point C, and requiring the local cross-section to be 
transverse also to the jet space (thus defining a cross-section for the prolonged 
action on J”), then one can find a moving frame for the action of the group on the 
multispace with the desired property, that is, the resulting moving frame will be 
the standard continuous moving frame when restricted to jets, and the discrete one 
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when restricted to M°"+)), The overarching continuity of the multispace manifold 
guarantees that one is the continuous limit of the other as the points coalesce. 
We note several features of Olver’s multi-space for pointed curves: 


(i) the curve C is not parametrised and the points z; on the curve need not be 

labelled in order with respect to some parametrisation, 

(ii) coalescence can take place between any two of the z; on the curve, see Figure 
10, 

(iii) none of the z; are distinguished in the sense that one of them is a natural base 
point for a projection of the multispace to M, 

(iv) the pointed curve j,C is essentially a set of points with a contact condition 
at each point. 


* O-jet 
1-jet 
Z2=23 Z2=2; 
Zo Zo 
Zı=Z4 
Z4 Z5 
Z5 Before coalescing After coalescing 


FIGURE 10. Coalescence of points in a pointed curve does not 
require the points to be adjacent, in some sense, on the curve. 
Note that there needs to be a well defined tangent of C at the 
point of coalescence. 


In our construction of a higher dimensional lattice based multispace, only a 
version of property (iv) remains. 

Our next section describes an interpolation scheme which can be applied to our 
geometric construction. 


4.2. Multivariate interpolation. One of the main problems with multivariate 
interpolation is that the solution to the interpolation problem is not unique in 
general and it might not even exist. For example, if we fix the values of a function 
f(x,y) at the two points (1,2), (—1,1) and we want to find a polynomial in x and 
y of minimum order, such that it coincides with the function at those points, we 
can use f(x,y) = a+ be or f(x,y) = a + by, and there is no reason why we would 
chose one over the other. On the other hand, if we fix the value of the function 
at (1,2) and (1,—1), then the first choice is not appropriate unless the function 
has the same value at both points, while the second one works. Thus, the choice 
of interpolating polynomial might depend on the data, it might not be unique, or 
even exist, and sometimes there is no reason to favor one choice over a different 
one. In the p = 1 interpolation case none of these problems exist. 

Thus, a main question in multivariate interpolation is: is there a family of polyno- 
mials which can interpolate the values of a given function and which has properties 
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like generality, minimal degree, uniqueness and having a well-defined Hermite poly- 
nomial (one for which not only the values of the polynomial, but also its derivatives 
at the points coincide with those of the function) as points coalesce? These are the 
properties we will need if we want to use them to define smooth coordinates in our 
multi-space. This question was answered by C. de Boor and A. Ron in [5], [6] and 
[7]; we describe below their solution to the interpolation problem as it applies to 
our particular case. 

Interpolation Notation Let II be the set of p-variate polynomials, and II’ its 
algebraic dual. Let A be a subspace of II’. We will denote by Il, the subset of 
polynomials of order less than or equal to k. 


Definition 4.3. We say that P C II is correct for A if for any continuous linear 
functional F on A there exist a unique q € P such that 


F(A) = Ala) 
for any A € A. We also say P interpolates A. 


The space II’ can be identified with formal power series at the origin o, using 
the bilinear form 


k (0) k oO 
(12) (af) = S32 oy 


k! 
See [5], [6] for more details. 


Example 4.4. Lagrange Interpolation. If A is spanned by point-evaluations, A = 
{Ao}oco with Ag(p) = p(0), 0 € RP, finding P correct for A solves a Lagrange 
interpolation problem. Indeed, if © represents a finite number of points in the 
parameter space D C R?, and og is evaluation at an element of O, then one can 
check that the power series representing g is the Taylor expansion of e?* (see [5], 
[6]). An analytic function F defines a continuous linear functional on A via (12) 
and F(\9) = (F,e**) = F(0). Thus, P is correct if for any F there exists p € P 
such that 


F (Xo) = F(8) = ro(p) = p(0) 
for all 6 € ©, which is the definition of Lagrange interpolation, see Example 4.9. 


Example 4.5. Hermite Interpolation. If we choose instead the set 


A = {Aqo}ocegevcn 


where 
q,0(P) = ro (a(V)(p)) = a(V) (p) (0) 

where V is the gradient vector and q € V C II is a properly chosen polynomial with 
coefficients reflecting the differential data we need to match (for example, q(x, y) = 
x? if we wish to match the second derivative in x), then the associated power 
series representing A,,9 is q(x)e®™ (see[5], [6]). Therefore, finding P correct for A is 
equivalent to finding a family of polynomials P such that for any linear functional 
on A represented by an analytic function F’, there exists a unique polynomial p € P 
with the property 


F(Aqo) = (F, (x)=) = a(V)(F)(9) = Aq,o(P) = a(V)(p) (9) 


for all 0 € © and all q € V. The different choices of V allow us to find unique poly- 
nomials that coincide with F and different choices of derivatives, or combinations 
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of derivatives, on ©. This is the solution to the Hermite interpolation problem, see 
Example 4.10.. 


Definition 4.6. Given a formal power series at x = 0, call it f, we denote by fy 
the homogeneous term in the power expansion of f of lowest order. 


For example, if 0,x € R?, efx =i. 


Definition 4.7. Given a finite subset A C II’, we can identify each of its elements 
with formal power series at the origin using (12), and we can consider A to be the 
vector space spanned by A as represented by these series. Define the vector space 


Ay := span{ f}, f € A}. 


For example, if we consider the planar case, and © = {(1,—1),(0,2)}, then 
eb-D)* =lta—y+... and e@2)* =142y4.... Since eb-D-* — @(0.2)* = 
x — 3y+..., we have that A, is the linear space spanned by the polynomials 
pi(a,y) =1= ey and po(x,y) = x — 3y = (eft. 1) — en yas 

From now on, if © is a data set, we say that A, is continuous on O if whenever 
the data ©’ is close to © in the standard product topology, then A is close to Ay 
in the standard topology of polynomial spaces. 

The following theorem is a compilation of results found in [7]. Since our con- 
structions are lattice based, we assume that A is defined as in the Lagrangian 
interpolation associated to a lattice of points © C R?, in which case we also write 
A as A®. We omit the superscript © where the dependence of A on @ is clear. 


Theorem 4.8. (|7]) The space A, has the following properties: 


ell defined. For any finite ©, the assignment A > exists, is unique 
1) Well defined. Fy y finite O, th ig t A > Ay exists, i ) 
and A, is correct for A. 
(2) Continuity. Recall that Il, is the set of p-variate polynomials of degree 
less or equal to k. If I, C Ay C Ig41, then the assignment A > A, is 
continuous with respect to O. 
oalescence = > Osculation. at is, the Lagrange interpolation be- 
3) Coal Osculati That is, the Lagrange int lation b 
comes the Hermite interpolation under coalescence, provided the coalescence 
is well controlled, so that data points coalesce along embedded curves. 
(4) Ay is closed under differentiation and it is spanned by a homogeneous basis. 
inimal degree. as minimal degree. 
5) Minimal deg A, h inimal deg 
(6) Monotonoticy. If © c O', then AP iG AF. 
(7) Cartesian product => tensor product. If © and ©' are two data 
sets, then 


(13) A = 19 @ AS’, 


(8) Constructible. The space A, can be constructed in finitely many arith- 
metic steps. 


Notice that once an interpolating family is chosen, the actual interpolation prob- 
lem reduces to solving a linear system of equations. Indeed, one would choose a 
linear combination of a basis generating A, and write a linear system for the coef- 
ficients using the values of the function we wish to interpolate on the interpolating 
data. The solution of the linear system will define the proper combination of the 
basis and hence the interpolating polynomial for the function. 
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It is essential that our construction of multispace ensures that the interpolation 
problem satisfies Properties (2) and (3) of the above Theorem. To quote de Boor 
and Ron [7], concerning Property (2) in the above Theorem, (note that A, is denoted 
as IIo in the original [7]) 

“I£ © c R? consists of three points, then one would choose Ie C 
II, (as our scheme does) but if one of the three points approaches 
some point between the two other points, this choice has to change 
in the limit and hence cannot change continuously. As it turns out, 
our scheme is continuous at every © for which II, C Hoe C Hk4+1, 
for some k.” 

Next, we quote de Boor and Ron [7], concerning Property (3) in the above 
Theorem. 

“ _..If, eg, a point spirals in on another, then we cannot hope 
for osculation. But if, eg, one point approaches another along a 
straight line, then we are entitled to obtain, in the limit, a match 
at that point also of the directional derivative in the direction of 
that line.” 

These limitations on continuity and coalescence mean that in our construction 
of our multispace, we cannot be as free in our choice of generalisation of the one 
dimensional pointed curves used to construct Olver’s one dimensional multispace, 
as might seem possible. We return to this discussion in §4.3. 


Example 4.9. In the Lagrange interpolation case, assume © C RP is given, as be- 
fore, by (1,2) and (—1,1), and assume {Xo }oco are the associated point-evaluation 
functionals. Thus, we have two series generating A, namely 

ety = 1444 2y+o(|x]), e77 =1—2+4+y+4 olx]. 
A basis for the vector space generated by these two series are the series of e™t?¥ = 
1+a+2y + o(|x|) and e?+2¥ — e~*t+¥ = 22 + y + o(|x|), and so 


Ay = span{1, 2x + y}. 
If we choose as © the points (1,2) and (1,—1), then the generators of À are the 
same as before, but 
Ay = span{1, y} 
as expected. 
If we choose four points of the form (1,2), (1+., 2), (—1,1), (—1,1+ 6), then the 
four series generating A are 


fi = 14 (1te)at2yt5 ((14-6)2t2y)?+0(IxP), fo = 1-#+(1+6)yt+-5 (—a+(L+e)y)*+0([|x|?). 


and 
L(+ 2y)? + ox), fa=1—atyt 5 (ety)? + olx). 


Since (fi — fs), = ex and (fo — fa), = ey, if € #0, then 1,x and y will be three of 
four generators for Aj. A fourth will be given by 


(Fe -f)+h-f)- hth) =@+0e- Fy, 


£ J 
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and hence 


i 
Ay = span{1, z, y, (2 + e)a? — >“ y?} 
will generate the interpolating polynomials. 
Example 4.10. In the Hermite interpolation case, assume we would like to instead 
find interpolating polynomials that coincide with a function at (1,2) and (—1, 1), 
and, say, with its partial with respect to x at (1,2) and with respect to y at (—1,1). 
In this case, the polynomials q; generating the Hermite data are, at (1,2), qı (£, y) = 
x, and at (—1,1) qo(x,y) =y. One can see (see [5], [6]) that rg, (1,2) is represented 
by the analytic function f(x) = xet?” , while Aq. (-1,1) is defined by g(x) = ye**¥. 
Since A = span{A 1,2); A(=1,1) Àq ,(1,2) Aqo,(—1,1) } ‘A has four generators, namely 
fi = xet? = g +g? + 2yrto(|x|’), fa = yet =y—ayty? + olll’). 


and 


E 1 ae 1 
fa = ®t? = l4r+2y+5 (2+2) +o(lx]?), fa =e 84 = 1-at+yt5(—aty)’+o([xI”). 


Some simple and direct calculations show that (f3), = 1, (f1), = x, (f2), = y and 
(2fi + f2 — fs + fa)y = 22° — Sy”. Thus, 


1 
Ay = span{1, z, Y, 2x? zj au} 
The interpolating polynomials in the previous example converge to these as € + 0. 
We are now ready to define a lattice based multi-space in several variables. 


4.3. Multispaces in several variables. We define two related versions of multi- 
space, the first containing the jet bundle J’(M,R) and the second containing the 
jet bundle J‘(U, M) where U C R? is open. 

We first recall that a point in the jet bundle J’(M,R) is represented by a triple 
(x, f,U], where x € U C M, the set U is open, and f : U > R, is a C” function. 
We say that the triple |x, f,U], ~ [a’, f’,U’], if « = x’ and if, in some coordinate 
chart containing x, f and f’ have the same derivatives up to order r, ([26], page 
60.) The equivalence class [x, f,U], is known as the r-jet of f at x. If T,(f)(x) 
is the order r Taylor polynomial of the (sufficiently smooth) function f at x, then 
(x, f,U], ~ [x,T,(f)(x), U]r, so we speak of T,( f)(x) as being the r-jet of f at a. 
Further, the coefficients of the rth order Taylor polynomials form local coordinates 
of the jet bundle J”(M,R). It is this construction that we generalise first. 

To construct our multi-space which both contains and generalises the jet space 
J°(M,R), we proceed as follows: 


1. We first define the kinds of lattices [ that we will take as the models of domains 
for a mixed discrete-continuous jet at a point in M. They will be sets of points in 
ZP with “directional multiplicities” or more precisely, “required contact conditions” 
attached. Our model lattices will come equipped with a base point. We show 
further that these models have the properties required for the de Boor and Ron 
interpolation of functions on them to be smooth, both as their image in M is 
varied and under coalescence. 

2. Next, for a model lattice T C ZP C RP, we let U C RP be an open set, 


diffeomorphic to the unit disc in R?, containing I. Let ¢ be a diffeomorphic map of 
U into M, and f : ¢(U) > R a function. Our multi-jet will then be an equivalence 
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class of quadruples |I, f,¢, U], where [[, f,¢,U] ~ [I’, f’, ¢’, U'], if the base points 
of ¢(T) and ¢(I”) agree; if o(T) = ¢/(T”) as sets; if whenever uz = o(x7) = G(x’ y) 
the required contact conditions on x; and z’ 3 are the same; and if the contact 
conditions induced on ¢(I) by those on T are all zero when evaluated on f — f’, 
see Figure 11. 

3. Finally, the multi-space fibre over x will be defined as the union of all equivalence 
classes of multi-space jets with base point x. Our coordinates on the fibre are 
those which assign to each |I, f,¢,U], both the coefficients of the image of T and 
the coefficients of the de Boor and Ron interpolant polynomial. In this way, we 
have the usual bundle topology on our multi-space which relates naturally to both 
local coordinates on M and to the coordinates on the fibre over x. Our multi 


space contains the jet bundle J‘(M,R) for each £, as an embedded submanifold, 
+ 
P ) 
(the number of derivative terms up to order £ on p-space), and the interpolation is 


given by the @th order Taylor polynomial. 


specifically as multi-jets where the lattice is a single point with multiplicity ( 


k a 
Xy -> XJ 
coalesce =(1,0 
a = dr w=(1,0) 
w 
w XJ xe Xj 2 
IR2 
M 
Ø Ø(U) 


FIGURE 11. After coalescence, the zeroth order contact condition 
at ug is replaced by D(p(f))lu,(v) = D(f)lu,(v) where v = 
D($)|e; (w). 


The second multispace we will define, containing J (Q, M) where Q is an open set 
of R?, is related to the first, by considering the function f, in the above construction, 
to be each of the coordinate functions on M, evaluated on the image of U. 


4.3.1. Basic definitions. As before, from now on we will assume that our lattices 
are connected. 


Definition 4.11. We say the lattice IT has an ¢-corner distribution, or a corner 
distribution of length £, if it has the following inductive description: 

If p = 1, the lattice is a connected lattice with + 1 vertices. Notice that the 
number £ refers to the degree of the derivative one gets when all points coalesce 
into one point, not to the number of points in the sublattice. 
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For any p the lattice is a connected lattice containing + 1, (p — 1)-dimensional 
disjoint corner lattices of lengths 0, 1, ..., 4. Figure 12 shows four corner lattices 
for p = 2 of lengths 4, 3, 2 and 2 (clockwise from the first quadrant). Figure 13 
shows a forward d = 3 corner distribution. We will consider corner lattices with a 
distinguished point xo. 


FIGURE 12. Corner lattices will allow the definition of finite dif- 
ferences needed to approximate a Taylor polynomial. 


FIGURE 13. A forward p = 3 corner lattice of length £ = 3. 


Corner distributions contain exactly enough points to define interpolating poly- 
nomials that will converge to Taylor polynomials upon coalescing. For example, if 
p = 2, and £ = 2, and uo is the base point, a possible interpolating polynomial 
will have coefficients which are a linear combination of the terms 


1, Ax f (uo,0); Ay f (u0,0), AZ. f (uo,0); AZ f(uo,0), Ay Az f(uo,0) 


where Agsluo o is the operator Az|uoo(f) = f(u1.0) — f(uoo), and similarly with 
Ay|uo,o- To be able to uniquely determine an interpolating polynomial with those 
coefficients we will need to use all the vertices in a corner distribution like the one 
in the first quadrant of figure 12, with length £ = 2 instead of 3. Different corner 
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FIGURE 14. Corner lattices under repeated coalescence of hyper- 
planes. An arrow indicates a zeroth and a first order contact condi- 
tion is required at that point in the interpolation, an arc indicates 
a zeroth, first and second order contact condition, a plane indicates 
all zeroth, first and second order contact conditions in the plane 
are required. The squared points are the base points. 


lattices will produce different types of interpolating polynomials, using forward, 
backwards or other types of differences. 

From now on we will assume that all lattices [ have a corner distribution of 
length £ with base point 29. When several lattices are involved we will denote the 
base point of [ by xr and the base point of (T) by ur = (zr). 


Lemma 4.12. Assume our data O is given by evaluating a function on the points 
of a corner lattice I of length £. Then 


A, =I 


Proof. Consider the monomials corresponding to polynomials of order 4. That 
is all monomials of the form z7 where |I| < £ is a multi-index I = (i1,...,ip) 
and a! = a'...a,?. We want to show that these monomials are generators for 
A. Thus, we want to show that they generate f}, where f is any possible linear 
combination of e”*, 6 € ©. Therefore, it suffices to show that the coefficients of the 
monomials x’, |I| < £, in the Taylor expansions of e”? form an invertible matrix. 
Notice that we have the same number of monomials as points and hence the matrix 
is square. 

Next, notice that the coefficients of these Taylor expansions are multiples of the 
monomials themselves evaluated at the point 6 (since we are simply substituting 
x; by Oix; in the expansion). Therefore, the matrix of coefficients is given by a 
multiple of the matrix with rows (6/), 0 € ©, where 07 has the different monomials 
in some prescribed order. This means that if the matrix were not to be invertible, 
we would have a combination 57, a 10! = 0 for a; € R which will be valid for all 
0 € O. That is, the points in © lie in an algebraic variety of order £. 
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But this is not possible: our points lie on £ + 1 distinct hyperplanes, and £ 
of them contain enough points to make the hyperplane unique (only one of them 
contains a single point and does not determine it). That means the polynomial 
must factor through the £ linear equations that define the hyperplanes and must 
have at least order £. But the extra single point left does not belong to any of 
the hyperplanes, and hence to ensure the point also lies in the variety we will need 
to use a polynomial of order at least £ + 1. Therefore the matrix of coefficients is 
invertible and 

Ay = II,. 


This lemma ensures that property (2) in Theorem 4.8 is satisfied when we use 
corner lattices. Next, we consider coalescence of the point in the lattice, leading 
to a change in the lattice, an increase in the contact order and to our Lagrange 
interpolation becoming a Hermite interpolation. We restrict the coalesce to be along 
coordinate hyperplanes in the model lattices; we call these kinds of coalescence 
hyperplane coalescence. See Figure 14. Forbidden coalescences are also illustrated 
in Figure 15. We note that hyperplane coalescence maintains the basepoint of the 
lattice, although not in general the contact condition there. Coalescence means, in 
effect, that we consider some lattice points to not be distinct. 


H-A 
E- £6 


FIGURE 15. Examples of forbidden coalescence. We restrict coa- 


lescence to being along coordinate hyperplanes, which maintains a 
coordinate structure. 


Since we want the multispace to be closed under hyperplane coalescence, we 
consider coalesced model lattices to come equipped with certain required contact 
conditions specified at particular points of the lattice. Suppose two points, uo and 
u1 = uo + hv coalesce as h > 0. Then the interpolation goes from matching the 
values of the function f at ug and u1, to matching the values f(uo) and D(f)|u,(v). 
See Figure 11. 
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Repeated coalescence leads to higher order Hermite interpolation problems. If 
you begin with a corner lattice and coalesce along hyperplanes repeatedly, you will 
arrive at a single point at which which the interpolation for a function f is simply 
the Taylor polynomial for f, with the order of the Taylor polynomial being the 
length of the corner lattice. In this way, the jet bundle J“(M,R) is a subset of our 
multispace. See Figure 14. 

From now on we will abuse the notation and denote the (coalesced) model lattice 
also as I and the contact conditions on it by C(T). Under the map 6: U DT > M, 
o(T) inherits contact conditions which we will denote by C(#(T)). 

The result of Lemma 4.12 remains unchanged under coalescence, as shown in 
the next lemma. 


Lemma 4.13. Assume IT is a hyperplane coalesced corner lattice of length £. Then 
Ay =I. 


Proof. We can show this by induction. Assume that 6; = 0; + he; approaches 6; 
along the e; direction, that is, as h > 0, for j € J, (that is for those data points 
in a hyperplane). Except for the coalescing of a corner, we would have more than 
one point coalescing since we are using limits of hyperplanes. 

Once more we want to prove that the coefficients of the monomials x7 in the 
expansion of exp(@- x), for any 0 Æ 6; and a; exp(6;-x) for any j € J define an 
invertible matrix. As in the proof of Lemma 4.12, the coefficient of x7 in exp(0;- x) 
is given by 0} /I!, and one can directly check that the coefficient of x! in x; exp(0;-x) 
is given by 

1 ; 

nO “ki 
if I = (kı,..., kp) and I — e; = (kı1,..., ki —1,...,kp), for j € J. Assume that 
the matrix formed by these coefficients is not invertible. It means that there is a 
polynomial of degree £ of the form 


(14) 5 arx! 


[T|<é 


that vanishes on all 0 € © except for ;. As we saw in the proof of Lemma 4.12, 
we would need an ¢-order polynomial to describe a polynomial vanishing on those 
points, given by the product of £ linear equations representing the £ hyperplanes 
(one of them will not be unique if it corresponds to a corner, since it is determined 
only by one point, but we can just make any choice). Now, if the lemma were not 
true, we would additionally need to satisfy the relation 


XO 610; “ki =0 


[T|<é 


where I = (ki,..., kp). This is simply the derivative of the polynomial (14) with 
respect to x;. Thus, the polynomial (14) needs to have order of contact 2 at the 
points on the hyperplane that coalesced, Typically one would need a second order 
polynomial along the hyperplane to achieve that, which would increase the degree 
of (14) to 2+ 1. Higher orders of contact would result from a higher number of 
coalescing hyperplanes and the order of the polynomial would generically increase 
accordingly, proving the Lemma. 
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coordinate 
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FIGURE 16. The data for a multispace element. The Figure shows 
functions f and f’ which agree on the contact conditions at the 
lattice points, ¢([). We have [T, f,¢,U] ~ I, f’, ¢, U]. 


The definition of elements of our multispace as equivalence classes of functions 
which agree on images of lattices, is illustrated in Figure 16. 


Definition 4.14. We define the multispace (M, R) as the set of equivalence 
classes of quadruples |L, f,¢,U], where (T,C(T)) is a (possibly hyperplane coa- 
lesced) corner lattice of length 4; U is an open set of R? for some p < dim(M), 
diffeomorphic to the unit disc and containing T; the map ¢@: U > M is an embed- 
ding of U into a single coordinate chart of M, and the map f : ¢(U) > R is smooth 
of order £. We say two quadruples are equivalent, [[, f,¢,U] ~ T”, f’,@’,U"], if 
the base points of d(T), @(I’) agree, that is ¢(ar) = ¢'(ar-), if 1) = ¢'(T'), 
C(o(T)) = C(¢'(T’)) and the contact conditions evaluate to zero on f — f’, that is, 
both f and f’ satisfy the same contact conditions on each vertex of (I). (We note 
the contact conditions are linear.) 


4.3.2. The Main Theorem: multispaces are manifolds which contain the jet bundle. 
Denote by Li(M , R) the subset of (M, R) given by non-coalesced lattices, that is 
a lattice with zero contact condition in every vertex. Denote also by Jf(M , R) the 
space of regular jets of maps from p-dimensional submanifolds of M to R. The 
main purpose of this section is to prove the following theorem. 


Theorem 4.15. Let M be a manifold of dimension m. Then, there exists a topol- 


ogy and a differential structure that makes (M, R)P into a smooth manifold of 
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dimension (m + Herr’), with the jet space J5(M, R) as a smooth submanifold and 


with LE(M, R) as an open submanifold. 


Before we start proving the theorem we recall that the main difficulty in defin- 
ing interpolating polynomials is determining the family of polynomials with which 
we choose to interpolate. Once this is determined, the actual interpolating coeff- 
cients are simply given by the solution of a linear system of equations defined by 
the equality conditions we need to satisfy at the chosen points. Therefore, they 
will change smoothly insofar as the linear system (and hence the family of inter- 
polating polynomials) changes continuously. Lemmas 4.12 and 4.13 together with 
Theorem 4.8 show that we are indeed in the smooth regime of the de Boor and Ron 
interpolation method. 


Proof. First of all, let us show that (M, R)? is a bundle over the manifold M. 
Let us call zp the base point of IT. The fibre over u € M is the set of equivalence 
classes, |I, f, ¢,U] where u satisfies (ap) = u. Define 7 : (M, R) — M to be the 
projection map 


(15) n(ID, f, , U]) = ¢(zr) =u. 


Let us first restrict to L6(M,R). In order to have a well defined coordinate 
labelling on the lattice, we need to describe an ordering on the vertices. The 
first point will always be the base point. Although any ordering will do, we can 
order them by induction: if the lattice is one dimensional we move from smallest 
coordinate to largest coordinate in R. If the lattice is two dimensional, then we 
order the hyperplanes from smallest length to largest length, and then order in each 
hyperplane as in the one dimensional case. In addition to its place in the ordering 
of the lattice points, each point in the lattice will have a coordinate in M. Since the 
image of the lattice lies in a single coordinate chart, these coordinates are consistent 
across the lattice, once the chart is designated. Thus the dimension of the set of 
embedded lattices is m(Pr. This set can be considered as a submanifold of MY, 


where N = Cc. with the subspace topology, and evidently has coordinates in 
terms of the coordinates on M. 

If the lattices are coalesced, the coalesced vertices will have repeated coordinates, 
as many as the multiplicity requires. The set of coalesced lattices will have measure 
zero since it is the level set of some functions depending on coordinates. 

If we now consider the relevant polynomial interpolation py of a (sufficiently 
smooth) function f on the embedded image of the corner lattice, we find that it can 
be described by oe coefficients. Since, by construction, [L, f,¢,U] ~ [, pf, $, U], 
we may add these coefficients to the coordinates of the equivalence class. That 
these coefficients can be used as coordinates follows from knowing that we are 
in the regime where the de Boor and Ron interpolation is unique and depends 
smoothly on the data, and from the fact that the coordinates of interpolations on 
the coalesced lattices are well behaved limit points of those on the non coalesced 
lattices. 

A simple counting tells us that we have, in total, (m+1)(? T) coordinates needed 
to specify the equivalence class in each element of the fibre. 

The bundle (M, R) is taken to be the disjoint union of the fibres, and 7 as 
in (15) is its projection. We may take the usual bundle topology given by the 
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smoothness of the local trivializations defined by the coordinates 
V c (MR) => RC) x ROP) 

which turns the space into a manifold. The map 7 is clearly smooth since the 

basepoint is simply the first element in the list of lattice points. 

We have already shown that Jf(M,R) C (M, R), indeed, it is embedded as 
the submanifold whose first set of coordinates is the diagonal A in M^, where 
N= Pe: We note the standard jet space is J‘(M,R) = J£ (M,R). 

Notice that if we perturb slightly an uncoalesced lattice I, the lattice will re- 
main uncoalesced, and so the subspace of classes of the form [T, f,¢,U], with T an 
uncoalesced corner lattice, is open. 


Remark 4.16. Although we have required the image of a lattice in M to be 
within a single coordinate chart of M, this restriction is perhaps not vital. The 
generalisation requires, firstly, keeping track of which chart as well as the coordinate 
given by the chart, in the lattice part of the coordinates for the multispace element. 
Secondly, it requires the construction of interpolations which agree on intersections 
of coordinate charts. 


We now construct our second multispace, (Q, M J where 2 C R’, p < dim M, 
which represents local approximations to embedded parametrized p-submanifolds 


in M. 


Definition 4.17. We define the multispace (Q, M ) to be the set of equivalence 
classes of triples [[,¢,U] where T c U C Q c R? is a (possibly coalesced) lattice 
of length 4; U is an open set of Q C R”, diffeomorphic to the unit disc, where 
@:U — M is smooth of order £ and where ¢(U) is contained in a single coordinate 
chart of M. We say two triples are equivalent, [[,¢,U] ~ [I’,¢’,U’] if the base 
points of I, I” agree; if d(T) = ¢(T’); if C(¢@(L)) = C(¢'(T")); and if d(U) and 
¢'(U") have the same order of contact (as submanifolds) at (T) as indicated by 
C(¢(T)). 

Theorem 4.18. Let M be a manifold of dimension m. There exists a topology and 
a differential structure that makes (Q, mM)? into a smooth manifold of dimension 
2m(? 2) : 

Proof. The proof of this theorem is almost identical to that of theorem 4.15. In 
local coordinates in M, we can write ¢ = (¢1,..-.,¢m). We would then apply the 
process in theorem 4.15 to produce m(? Bo coordinates that determine the lattice, 


plus m(? 9 coordinates determining the interpolating polynomial for øk with data 
I, for each k = 1,...,m. The remainder of the proof is identical to that of theorem 
4.15. 


The main purpose of our multispace construction is to show that a frame on a 
multispace is simultaneously a frame on the jet bundle and a frame on the set of all 
local lattice based discretisations. We now proceed to discuss how the group action 
on M induces a group action on multispace. We then show that a moving frame on 
multispace is simultaneously a smooth and a discrete frame, with the smooth frame 
being the limit of the discrete, and that the discrete Maurer Cartan invariants and 
their syzygies coalesce to the smooth ones. 
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4.3.3. The action of a group on (M, R). Let G be a group acting on M x R, 
GxMxR-MxR. 


Recall the equivalence classes in (M, R)? have the form |I, f,¢,U] where rT CU c 

R?, ¢ : U > M and f: d(U) > R. Then for each g € G there is an induced action 

on the map (¢, f) : U x 6(U) + M xR given by (g: (¢, f)) (x, z) = g: (0(2), F(2))- 
The action of G on (M, R) is given by 


g: IP, f$, U] = (0; (9: ($, F))2; (9° (6 f))a, U] 


We note that (g-(¢, f))1 may be an embedding only for g in a neighbourhood of 
the identity. In this case we would have a local group action as in [52]. Since the 
action of G on M x R preserves the order of contact, this action is independent of 
the representative of the class and is thus well-defined on (M, RM, 

If, in a particular application G acts only on M, one can extend to an action on 
M xR by taking g- (z,t) = (g-z,t), that is the identity action on the R coordinate. 

If [ is a single point so that [T, f,¢,U] is an element of JE(M, R), the induced 
action is the standard prolonged action, that is, as induced by the chain rule, while 
if the lattice [ is uncoalesced and the group does not act on the parameter R, the 
action is the diagonal action on MN where N = (ek 


4.3.4. The action of a group on (Q, muy. Let G be a group acting on M 
GxM—>M 


and let [[, 6, U] € (Q, M$ be an element of our multispace. Since ġ : U C Q > M, 
we may define g-¢: U > M by (g- $)(x) = g - (¢(a)). Define the action of G on 
M$? to be 
g IP, ¢,U]=[0,9- $, U]. 

Again, since the action of G on U x M preserves the order of contact, this action is 
independent of the representative of the class and is thus well-defined on (Q, M ). 
Further, the action restricted to Bpltve M) c (Q,M ) is the standard prolonged 
action induced by the chain rule. 

Notice that one could consider more general actions on M x U, but we will omit 
it here to avoid further complications. 


4.3.5. Moving frames on (M, R), (Q, mM). We are now in a familiar situation: 
we have a smooth manifold (M, RM, or (Q, M ) and the action of a Lie group 
G on it. Thus, we can investigate the use of the standard moving frames method 
developed in [17] to establish the existence of a moving frame for the multispace. 
Assume the action of the group G on (M, Ry, or (Q, My), is free and regular, 
so that the existence of a local moving frame is guaranteed. Let us choose a point 
in L € (M, R)? or (Q, mu)? and let S be a section transverse to the orbit of G 
through L. Using the standard moving frame method, we would get a local moving 
frame p, defined for all Ê € (M,R)$”, or (Q, my? in some neighbourhood of L, as 


the group element such that p- L € S. That is, p is an equivariant continuous map 
p:UC(M,R)O >G (resp. p :U C (Q, M) > G) 


where U is an open neighborhood of L. 
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Remark 4.19. We note that there are results detailing conditions under which 
an action on a jet bundle will become free and regular for a sufficiently large pro- 
longation, that is, by considering sufficiently high order derivative terms [17]. A 
discussion of the related results for a product action, under a sufficiently large num- 
ber of products, is given by Boutin [8]. We conjecture that similar results will hold 
for actions on multispace. 


Before proving our results, we give a simple example. The first example refers 
to the multispace, (R?, R)®. 


Example 4.20. We consider the two dimensional group G = Rt x R acting on 
MxR=R?xR as (\,€)-(z, y, u(x, y)) = (x,y, AU+e). We take ¢ to be the identity 
map for simplicity. At the corner lattice T = {(2o, yo), (£o + h, yo), (£0, yo + &)} 
with (xo, Yo) the base point, the multi jet coordinates are the lattice coordinates, and 
the coefficients of the linear interpolant to some function u on these three points. 
The interpolant is 


plu)(x,y) = A+B.(x—20)+C(y— yo) 
= sons Meth ay 
u(%o, Yo + n — u(zo, Yo) E 


Thus the coordinates coming from the interpolant are 
u(xo + h, yo) — u(zo, Yo) Uu(xo, yo + k) — u(zo, Yo 
(A, B,C) xs u(2o, Yo); l ) ( ) ( ) ( ) “ 
h k 
We see that the coefficients are functions of u at the lattice points and so the induced 


group action on these coordinates is the natural action on functions. We thus have 


(à, €) (A, B,C) = (AA + €, AB, AC). 


Theorem 4.21. Assume p is a local moving frame for the action of G on (M, R) 


(resp. (Q, Myf) determined by a section transverse to an orbit of G. Assume that 


the section is also transverse to the orbit through a point L € JMO (M, R) (resp. 
J°(Q,M)), that is, the domain of the multispace frame includes L. 

Denote by p(Q) the multispace frame at Q, and by p(L) the smooth moving frame 
which is obtained by the restriction of the normalisation equations for p to the jet 
space, evaluated at L. If L and Q have the same base point, then as Q coalesces to 
L 


d 


p(Q) > p(L). 


Proof. We prove the result for (M, R), the other case being similar. Notice 
that p(L) is the standard moving frame on jet spaces obtained through a section 
transverse to the prolonged orbits, which is the multispace section restricted to the 
jet bundle. We note that JÉ (M,R) is a submanifold of (M, RM, invariant under 
the action of the group, so that the orbit of G through a point L € JEM, R) is 
equal to the prolonged orbit. Furthermore, if a section is transverse to the orbit 
within (M, R), its intersection with J{(M,R) will also be transverse to the orbit 
within the submanifold I (M,R). 
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The proof is now immediate from the fact that the moving frame on (M, R) 


is a smooth map. 


In Appendix A, we will show a different convergence theorem for families of 
discrete frames. 


Example 4.22. Continuing with our previous example 4.20, a moving frame p is 
defined by (A,€) - (A,B,C) = (0,1,*) where * will be the invariant, (A,€) - Clp. 
This yields (A, €)|p = ( hu(zo.yo) ) or in the standard 


h 
u(xo+h,yo)—u(zo,y0)? u(xo+h,yo)—u(xo,Yo) 
matrix representation of this group, 


h hu(xo, Yo) 
(16) p= u(zo + h, yo) — u(xo, yo) u(ro + h, yo) — u(xo, yo) 
0 1 


It is readily seen that as h + 0, we obtain the correct induced group action on Uz, 
and that the limiting frame 


1 u 
(17) p(T, U, Ue) = | Ue Ue 
0 1 


is obtained both from the limit of the normalisation equations, (A, €); (uz, u) = (1,0) 
as well as being the limit of the frame itself. In other words, both the discrete and 
the smooth cases are handled by the one calculation of the frame on multispace. 


4.4. The continuous limit of invariants and syzygies. We return to our dis- 
cussion of discrete moving frames in the previous section, in which we have a lattice 
variety Ly embedded in some manifold, and a discrete frame is a map p : Ly — GN 
where N is the number of points in the lattice. Suppose now that adjacent ver- 
tices in the lattice variety begin to coalesce. Under what conditions will the discrete 
frame converge to a smooth frame, p? Furthermore, under what kinds of conditions 
will we have 


(18) Pnt1Pn ~ Plz +ev)p(z)* ~ exp (e (Dp) p*) 


where D is an invariant differential operator? And when will the discrete local 
syzygies converge to the local differential syzygies? 


Example 4.23. Let us go back to our running example 4.20. Setting 
pi = p(Ti, Yi, U(@s, Yi), U(Li + hi, yi), W(@e, Yi + ki)) and calculating pıpg', we have, 
using Equation (16), 


hı(u(xo + ho, yo) — u(z0, y0) hı(u(zo, yo) — u(z1, y1) 
pipo' = | Polul + hi, y1) — ulzi, y1)) ulti + hi, yr) — ulzi, y1) 
0 1 
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If we now set (1,91) = (£o +h, yo) and ho = hi = h (say), so that the second 
lattice is the shift of the first, by (h,0), then 


A |z 
Uz(zo + h, yo) — uz (xo, Yo) huz(zo, yo) 
= u(zo + h, yo) — u(Xo, Yo) u(ao + h, yo) — u(xo, yo) 
0 0 
_ bee =i 
+h 50 Ur (x0,Yo) 
0 0 
= —1 
= Pap |m) 


Alternatively, using the method we will apply in Example 4.27, we have setting 
u(x +h, y) = u(z, y) + hus(z, y) + O(h?), that 


(19) 
pıPo ` 
Ua (Xo +h, yo) — Ux(Zo, Yo) 7 T2 huz(x0, yo) F 72 
h+O(h h+O(h 
= u(ao + h, yo) — u(xo, yo) G u(zxo + h, yo) — u(xo, yo) a) 
0 1 
5 ugz(x£o + h, yo) — Ux(Xo, Yo) hux (xo, yo) 
~ exp [h u(zo + h, yo) — u(xo, yo) u(zo + h, yo) — u(xo, yo) 
0 0 


~ exp (hpzp"*) 

as above. We note that 0/Ox is an invariant operator since the independent vari- 
ables are invariant under the action, so that prp™! is invariant. Cases where the 
independent variables participate in the action require more care, as we indicate 
below. 


Our first theorem concerns the convergence of the discrete Maurer—Cartan in- 
variants to the smooth Maurer—Cartan invariants of the smooth frame. 

We consider the case where the discrete frame (pp) on the lattice variety Ly C M 
can be viewed as a multispace frame p, with 


ps = p([Fs,¢,U)), us = $(x3), 
where zy is the basepoint of Ty 
Assume the point Tixy is also part of Ty so that Tius = ¢(Tixz). If we have a 
path from 7;(uy) to wy in the multispace indicating their coalescence, we can use 
the multispace frame to differentiate 7; (p,7) es with respect to the path parameter 
at uz. 


Theorem 4.24. Let a multispace M with the embedded jet bundle J be given. Let 
a path in M be given, € > [F(e),¢,U] for0 <e <1. Let u(e) be the base point 
of o(T(e)). Assume that both paths and the coalescence (T (€)) > (To) lie in the 
domain of a multispace frame p. Set ple) = p((['(e), ¢,U]), that is, p evaluated at 
the point [T (€), p, U]. If v = d/de|-ou(e), then 


| POO] = (Devo) o* 


2 li os 
(20) rO)STo E 


e=0 
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where D(v) = X` v; /ðx; is the directional derivative. 

If the independent variables participate in the group action, then (D(v)p) p7} 
will not be invariant. We recall that a smooth frame on a jet bundle yields a 
canonical, maximal set of invariant differential operators. Indeed, on a manifold 
with coordinates x, if g- x = X, then we define 


o 


A K 
g=p Ox ij g=p OX ; 


Rewriting the partial derivatives in (20) in terms of the invariantized deriviatives, 
by inverting equation (21), yields an invariant expression which will satisfy equation 
(18). Indeed, suppose that d/de = D(v) = $0, v,0/0x;. Inverting Equation (21) 
yields expressions of the form 0/0x; = X, AikDk, so that 


(21) D, 


O O = Daw =0@ 
$ k k 


i k 


where the last equality defines the vector V, and D(V) is the invariantised directional 
derivative. 

The theorem follows from standard calculations with smooth moving frames. 
The invariantisation of the right hand side of Equation (20) can be calculated using 
the cross section that determines the frame and the infinitesimal action, see [37] for 
a detailed discussion. 


Example 4.25 (Special Euclidean group action on curves in the plane). The group 
is SO(2) x R? with the standard linear action of translation and rotation of curves 
on the plane, specifically, 


Cee Genie ea) 


This is the inverse (or right) action and so the equivariance of the frame is p(g-z) = 
gpz). 

The multispace frame calculation is as follows. We take an order 2 interpolation 
as we wish to achieve a multispace approximation of the curvature. If we interpolate 
the curve (a(t), u(t)) at T = {to, to + hi, to + h2} with base point to, we get 


p(x(t)) = A(x) + B(z).(t— to) + 3C (2).(t — to)? 


(h? — h2)x(to) + h2x(to + hi) — hix(to + he) 
hyha(h2 — hı) 

PUZ m hı)z(to) cme hox(to + hı) + hya(to + h2) (t 

hyha(h2 — hı) 


x(to) 


(t — to) 


to)? 


DISCRETE MOVING FRAMES Al 


and similarly for p(u(t)) = A(u) + B(u).(t — to) + $C(u).(t — to)?. The induced 
action on the coefficients is that induced on (a(t), u(t)) so for example 
(ht — h3)g-u(to) + hog: u(to + hi) — hig - ulto + h2) 

hyha(h2 — hı) 


g: B(u) 


1 
h2 
neea 


+h3(cos O(x(to + hy) = a) + sin 0(u(to + hy) = b))) 


h3)(cos 0(a(to) — a) + sin 0(u(to) — b)) 


—h? (cos O(x(to + h2) = a) + sin O(u(to + h2) = b))] 


The normalisation equations g- A(x) = 0, g- A(u) = 0 and g- B(u) = 0 yield the 
frame at T to be 


(hî — h3)u(to) + hõu(to + hi) — hfZu(to + h2) 


BEM aN UM) BONE Gore ae ts) © ea io tn) = hay Ete) 


In the limit as hg > hı, we have 


Ut (to hy) 2u(to hy) 2u(to) 
zlto hı) 2x (to hı) 2x(to) 


tan > 


and then finally as hy + 0, we have tan + u;/x, as expected, indeed, yielding the 
smooth frame as determined by the limit of the normalisation equations. 

If we calculate the invariant matrix p([I', u, ¢,U])~'((L, u, ¢, U]), recalling we 
have a right action (we are taking ¢ to be the identity map), we obtain the dis- 
crete multispace Maurer-Cartan invariants for this frame. Further, p- B(x) > 


(a? + uz)? so that we may treat p- B(x) as the multispace approximation to the 
infinitesimal arc length. 

In the above example, we used an invariant parameter t to describe the curve 
(a(t), u(t)). If instead we parametrise the curve as (x, u(x)), so that the parameter 
participates in the group action and the operator 0/Ox is not invariant, then greater 
care is required. For example, the group action on the interpolation curve, written 
as p(u)(x) = A(x, u) + B(a, u)(a — zo) + C(x, u)(x — x0)", looks like g - p(u)(x) = 
A(T, Ù) + B(T, U) (£ — T0) + C(&,u)(x — T0). Solving for the multispace frame in 
this case seems nontrivial. Such examples will be examined elsewhere. 


Now assume we have four lattice points, uy, huz, huz and JoTiuy = 1 Tou,, 
and that we have two paths connecting uy with 72Tjuy = TiTouz via each of 
the Tiuz and J2u, respectively. If we can associate the discrete frame with a 
multispace frame and differentiate the local syzygy associated with the discrete 
frame with respect to the path parameters at uz, we obtain the differential syzygy 
associated with the multispace frame at uz. Indeed, let a discrete frame (pJ) at 
the points uy, Tuy, Touz and 7, J2u,y be associated with the multispace frame p so 
that pz = p(zz) for some zy = [T 7, ġ, U] in the relevant multispace M, uy = ¢(a,). 
Then, the local syzygy is 


(22) T(K2)Kı = To(K1) Ko, Ky = p(Tizs)p(zz)*, Ka = p(Tazs)p(zs)*, 


where Tizy = [Til s,¢,U] and TiT y is the shift of Ty by Ti, with base point Tix, 
and Tjuy = $(7Tizz7) = O(@s+4e,). Here we are assuming that the shifts of Ty 
remain within U. Now let z‘,(e;) be paths in the multispace lying within U, with 
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25(0) = 27, 21) = Tizs, 25(e) = Ei (a), 9, U] and uf (ei) = o(x5(e)), v5 (a) the 
base point of T% (e;). We denote 


d d 
es de, wattle); = deg nots (@)s 
and 
Kes E| ola), Kw =] lee) le) 
deyle=o’ “7 ? dez leo! 7 


Theorem 4.26. After differentiating twice, once each with respect to the path 
parameters €, and ez and coalescencing the lattice of the multispace point zz to its 
associated jet point, and under the conditions just stated, the local syzygy, equation 
(22), becomes the continuous basic syzygy associated to p(zz) 


(23) D(w)Ky = D(v)Kw + [D(w), D(v)]pp~* + [Ky, Kw] 


Proof. The core of the proof is standard. We assume a matrix representation of 
the frame, and note that Taylor’s Theorem is valid. 

To ease the notation, we set €&1 = h and eg = k, and simplify to where v = 
d/dh|n=0¢(#o + he1) and w = d/dk|k=0¢(x£0 + ke2) in local coordinates, so that we 
evaluate our frame at the multispace elements with lattice basepoints at uy = (zo), 
Tuy = Uzi = O(@o + hey), Toug = usa = G(X + keg) and T2Tuy = usi = 
Us21 = (xo + he; + kez). We denote the partial derivative operators 0/0«; as Oj. 

Denoting Tip7 = pJ; we have 


Kja = p317", K352 = p32p7", 


and for sufficiently small h and k, there will exist, dropping the index J for clarity, 
matrices X; and Xə in the Lie algebra of G such that 


Kı = exp (hXı + O(h’)), Ko = exp (kX2 + O(k°)) . 
We have 
(24) 3p- p = Xi, 2p- p! = X2 
and 
ToX1 = Xı +kO2X1+O(h?),  TiX2 = Xə + kô2X2 + O(k?). 


Then 
Tok, = exp hToX1 T O(h?)) 
exp hXı T hkO2X 1 T O(k?, h?)) 
Ti Ko = exp kT, X2 =e O(k?)) 
= exp kXo T hkO, Xə T OR, h?)) 
Applying the Baker Campbell Haussdorff formula, [57], 


log (exp(X) exp(Y)) = X +Y + $[X, Y] + higher order brackets, 


II 


we have 
log (Tek, s Kə)) = hx, T hkð2Xı ZN kXə Ae 4 [hX + hkô2Xı, kXo] T O(k?, h?) 
= hXı +kXə + hk (32X1 + 4 [X1, X2]) + O(k?, h?), 
log (1K : Kı)) = kXo alg hki Xo hXı 5 [kX t hkO, Xo, hXı] a Ok, h?) 
= hXı +kXə + hk (31 Xə + 4 [X2, X1]) + O(k?, h?) 
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Equating the two formulae by imposing the local syzygy, differentiating with respect 
to both h and k and then sending h, k — 0, yields, after a slight rearrangement, 


O2X = OX m [X2, Xi] . 


(25) 


Finally, we need to rewrite Equation (25) in terms of the invariant differential 
operators. The formula given in Equation (21) shows that the partial derivatives 0; 
can be written as a linear sum of the invariant operators with invariant coefficients. 


—1 


We must then back substitute for the 0;, including rewriting the X; = O;p- p~* in 
terms of the Djp-p~'. The final result yields the extra terms in the case that the 
invariant operators do not commute. 


Example 4.27. We conclude the running example, 4.20. We set the points (£i, Yi), 
i = 1,2,3 to be (x1, y1) = (to +h, yo), (£2, y2) = (0, Yo + k) and 
(v3, y3) = (to + h, yo + k). We then calculate the four matrices, Kio = pipp’ = 
exp (Xio), Ks = p3p;' = exp(X31), Koo = popo’ = exp(X20) and K32 = 
p3p3" = exp (X32). Direct calculation gives, setting AF = F(ao+h, yo)—F (20, yo) 


for a function F € {u, Uz, Uy, Uzy} in the formulae to ease the notation, 


X10 


Aug Ux (0, Yo) 


_h Au Au ai O(h?) 
0 0 
Aiuy uy(To, Yo) 
iE Ayu Ayu + O(k?) 
0 0 


Xoq — hk Xai + O(R2, k?) 


AjuzAiuy + Ay Ugy Aru 


UyA1Ug + Uny Aru 


Xə — hk (Ayu)? 
0 
X19 — AKX32 + O(h?2, k?) 


AjtizAjuy + AjUgyAiu 


(Ayu)? 
0 


Ug Ai Uy + UsyAiu 


Xio — hk (Ayu)? 
0 


(Ayu)? 
0 


+ O(h?, k?) 


+ O(h?, k?) 


where this defines pet and Kiss The local syzygy is K31K19 = Ks2K20, and apply- 
ing the Baker Campbell Hausdorff formula to this yields 


The equation for the lower order terms simplifies to 


(26) 


1 
X19 + X31 + = 


2 


1 ae 
[X31, X10] = X32 + X20 + z X32; X20] + O(h?, k?). 


X32 — X31 = [X10, X20]. 


This last equation is straightforward to verify. Finally, taking the limit as h > 0, 
Equation (26) yields the differential syzygy for p evaluated on the jet bundle, 


o o 


-1 =ï —1 -1 
D Bp PvP = [oyp pap |, 


where recall p on the jet bundle is given in Equation (17). 
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Similar relationships exist when we take limits only in one of the variables, 
producing an evolution of discrete submanifolds. For example, if p = 2 and S' has 
coordinates (x,y), then if K} = pr+e.(z4)e(zz)~' and Ny = #£(Kj) Kj", then 
when we take limits in the calculation in the proof of the theorem as 7j4., > TJ 
we have 


d = 
ag EoD) = Notes — Kua NIK Ga) 


which is a mixed syzygy that often appears describing invariant evolutions of poly- 
gons in terms of coordinates in their moduli spaces, as in [41] and [46]. Among 
these evolutions one often finds completely integrable discretizations of well known 
completely integrable PDEs. These results are really key to some of the applications 
in our next section. 


5. APPLICATIONS 


5.1. Application to the design of a variational numerical scheme for a 
shallow water system. This example is motivated by the need for finite difference 
versions of variational shallow water problems which are invariant under the so- 
called particle relabelling symmetry. We consider the base space to have coordinates 
(a, b,t), where (a,b) is the fluid particle label at time t = 0,. The two dimensional 
dependent variable space is (x, y) = (x(a, b, t), y(a, b, t)), which is the position of the 
fluid particle at time t, so that (x(a, b, 0), y(a, 6,0)) = (a,b). The particle relabelling 
action is given by 


g-a=A(a,b), g:b= B(a,b), A,B, — ABa = 1 


together with g- x = x, g: y = y and g-t = t. It can be see that the particle 
relabelling group is the group of area preserving diffeomorphisms of the (a, b) plane 
(or at least the domain of interest in the (a,b) plane). Further, it is known that 
the invariants of this group action are A = ray» — XpYq and its derivatives under 
the invariant differential operators, 


O WO Yad ð at, 0 rað 


by Ada A Oe “Oy Ada A Ob. 


The aim is to design a multispace version of the Lagrangian for variational shal- 
low water problems, which have the form, 


(27) Liz, y] = Jiven tan — TpYa) da db dt. 


This family of Lagrangians are each invariant under translations in time, transla- 
tions in both a and b, rotations in the (a,b) plane, and more generally, the full 
particle relabelling group. Noether’s Theorem [48, 1, 52] then yields conservation 
of energy, linear and angular momenta, and potential vorticity. 

If we take the simplest corner lattice with base point (ao, bo, to) to be 


T = {(ao, bo, to), (a1, b1, to), (a2, b2, to), (as, b3, t1)} 
then the (linear) interpolation of x is given by 


x(a, b, t) ~ x(ao, bo, to) + M(x) (a — ag) + M(xy)(b — bo) + M (xz) (t — to) 
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where this defines the coefficients M (xg), and we have 


1/1 £(a9,b9,to) bo 1/1 æ x(a, bo, to) 
M(ta) =F 1 w(ai,b1,to) bi |, M(x) = 7 1 a, 2(a1,61, to) 
1 x(az, b2, to) bə 1 ag x(ag, b2, to) 
and 
1 ao bo x(a, bo, to) 
1 1 ay by x(a by to) 
M = 9 i 
(z+) (tı— to) A| 1 az b2 zx(az,b2,to) 
1 a3 b3 x(az,b3,t2) 
where 
1 ao bo 
A= 1 ay by 
1 a2 bo 


is the area, |(a@, — ao, bı — bo) A (a2 — ao, b2 — bo)|. The interpolant for y is similar, 
with y(a;, bi, ti) replacing x(a;, bi, ti) in the above formulae. 

We consider the Lie group SL(2) x R? acting on the (a,b) plane as the standard 
(right) equiaffine action, 


a g:a\_ ô -8 a— € 
wl ie ae ae) 
so that (€1,€2) € R? is the translation vector, and aô — By = 1, and that g: x = 2, 
g:y and g-t=t. This group is contained within the particle labelling symmetry 


group, and is just big enough to obtain the area invariant, which we do next. 
The induced action on the coefficients in the interpolants is given by, for example, 


1 1 x(ao, bo, to) Gg: bo 
g: M(za) = 1 x(aı, bı, to) g: by ; 
1 x(az, be, to) g: be 


noting that g- A = A, indeed, A is an invariant as is easily seen. 
We take the normalisation equations g - (ao, bo) = (0,0), g- M(#a) = 1, g- 
M(za) = 0 and g- M(ya) = 0. Then the multispace frame M(p)(T) is (€1, €2) = 


(ao, bo) and 
( 5 £)-( Be ae) ) 


T ee MA)  M(A) 


where M(A) = M(za)M(yb)— M(z)M (Ya). Evaluating the remaining coefficient 
on the frame, we obtain the invariant, 


I(M(y)) = g: Mv) 


M(A) 
1 Teg ao, bo, to) y(ao, bo, to) 
— A 1 a ay, b1, t1) ylar, 1, t1) 
1 x(a2,b2,t2) y(az, be, t2) 
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Calculating the continuum limit of M(A) we obtain tayp — £bYa, which is A, the 
area invariant, as expected. Further, the continuum list of the frame is 
Ta Th 


p> Ya Yb 
LaYb — LoYa ZaYb — XoYa 


This is the smooth frame obtained with the smooth limit of the normalisation 
equations, that is, {£a = 1,7, = 0, Ja = 0}. 
We observe that both M (x+) and M (y+) are invariant under the equiaffine action. 


PN 


FIGURE 17. Length one lattices stacked as a mesh, for a finite 
difference variational problem. Shown here is single time slice for 
the shallow water problem. 


Thus, we propose the multispace analogue of the Lagrangian (27) to be 


(28) M (Lix, y]) = XC L (x,y, M(z4), M(y), M(A)) A (tr — to) 
r 


where the sum is over all corner lattices stacked into a mesh, as in Figure 17. The 
factor A (tı — to) is the multispace approximation of the volume form, dadbdt, and 
is needed to obtain the correct continuum limits for the conservation laws for energy 
and the linear momenta. 

Finite difference Euler Lagrange equations and Noether’s conservation laws can 
be calculated in the standard way [32, 23]; the details and the results of this calcu- 
lation will be explored elsewhere. Here we simply demonstrate that the multispace 
Lagrangian, (28) is invariant under a discrete analogue of the particle labelling sym- 
metry. Indeed, looking at Figure 17, one can use a different element of SL(2) x R? 
on each corner lattice, provided that certain consistency conditions hold, specifi- 
cally, that if a vertex is in the intersection of two lattices, that their image under the 
two different group elements is the same. The Lagrangian (28) is clearly invariant 
under this discretisation of the particle relabelling group, the discretisation being 
a subgroup of (SL(2) x R2)* where R is the number of corner lattices on a time 
slice of the mesh. Finally, if one takes the collection of frames consisting of the 
multispace frame restricted to each vertex, then we have a discrete moving frame, 
using which the discrete Euler Lagrange equations and discrete conservation laws 
can be studied [53]. 
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5.2. Discretizations of completely integrable PDEs. The geometry of curves 
and surfaces have been linked to integrable systems repeatedly in the literature, 
see [24], [56], [44], for example. A drawback of the application of the results in 
this paper to finding completely integrable disccretizations of completely integrable 
systems is that one needs to choose a type of approximation (forward, backwards, 
linear or higher order, etc) a priori to find the limit. On the other hand, any 
two choices of discrete moving frames (be the one associated to a certain type of 
limit or any different one) will always be associated by a gauge transformation. 
This means that if one finds a discrete integrable system associated to any given 
choice of moving frame, one might be able to relate it to a different choice and 
perhaps link it to the continuous case. This was done in [46], where the authors 
found discrete integrable systems that were the discretization of Adler-Gel’fand- 
Dikii integrable evolutions, both of them linked to the projective geometry of curves 
and polygons. The authors of [46] also found a way to obtain two Hamiltonian 
structures associated to the discrete system through a reduction process, a process 
that was later extended to other semisimple homogeneous spaces in [45]. Different 
approaches were used, for example in [2] and in other works of these authors to 
construct completely integrable discretizations of integrable systems with the use 
of lattice models in Euclidean, projective and conformal geometry. Their approach 
is quite different from the one used in [46] and in this paper, in that they choose 
lattices with different geometric properties to achieve discretizations. The following 
are just some examples of the connection between continuous and discrete models, 
in the equi-centro affine and the Euclidean space. Both of these examples use mixed 
multispace discrete/continuous models, and although the first example only uses 
discrete coordinates to construct the multispace moving frame, the second example 
will make full use of both. 


5.2.1. Integrable discretizations of Boussinesg equations. In this example we make 
use of mixed differential/difference coordinates in the multispace. For reasons that 
will become clear later, we will also assume that the lattice variety has a monodromy 
in the discrete variable (a global property). That is, un+r(t) = m - un(t) for any 
n € Z, with T the period and some monodromy m € SL(3,R). This ensures that 
the invariants will be T-periodic in n. 

First we describe the situation when we are in continuous jets with two param- 
eters (x,t). It is well known (see for example [9]) that the Boussinesq equations 


1 2 
(qo)t + gi. + 304 =0, (qı): — 2g) = 0 


can be obtained as the evolution induced on equi-centro-affine curvatures by a 
certain evolution of curves. Let our manifold be R? with M = SL(3,R) acting 
linearly on it. Within M consider parametrized surfaces on (x,t). We will define 
the following cross-section: 


(29) p:u=e3, pu =e, pu” = de 


where the prime denotes d/dz and e; are the standard unit vectors in R3. Clearly 
d = det(u, u’, u”). This defines uniquely a right moving frame whose left companion 
is given by 
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Let us assume that 
(30) d = det(u”,u’,u) = 1, 


that is we will parametrize the surface so that the curves associated with t fixed are 
parametrized by the equi-centro-affine arc-length (these curves are in one-to-one 
correspondence with projective curves, see [9]. In that case the left «-Maurer— 
Cartan matrix associated to it is given by 


010 
Q= = |k 0 1 
ko 0 0 


where u” = kyu’ + kou. Next we will gauge this frame to a different left frame 
P =p 'g by the element 


1 00 
g=| 0 10 
—kı 0 1 
The resulting x-Maurer-Cartan matrix is given by 
0 1 0 
K=g'q¢ +g 'Qg=|0 0 1 
b a 0 


where a = kı and b = ky — ki. Gauging the system can be seen as changing the 
coordinates, the results can always be gauged back to the original setting. 
We will next consider the syzygy 


2 
(31) u — u” + zou = 0 


which describes a precisely chosen evolutionary equation for curves whose flow will 
be tracing our parametrized surface. With this condition, the left t-Maurer-Cartan 
matrix is easily seen to be given by 


—(wı +1/3a) wo 1 
N= OB = UL Wi 0 
U2 W2 ža 


for some entries v;, wi. The local basic syzygies (or the compatibility condition 
between x and t) are given by 

Ki = N; + [K, N] 
and they can be used to solve for N so that 


wo = 0, wi = 34 wa = b+ aa’, w=bt za, vg = —b' — =a 


We can further find two more syzygies given by 
2 2 
a, — 2b’ — a” =0, bi + 32 +" — 30% 
This system of equations is equivalent to the Boussinesq equation. Indeed, if a = q 
and b = 4q/, — qo, we have 


1 2 
(go) + ga + xan = 0, (qı): — 249 = 0 


which is the standard Boussinesq equation. 
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Assume now that we move in the multispace away from a continuous jet to a 
mixed discrete/continuous multispace submanifold, where x is now discrete and t 
is continuous. Let us choose lattices containing (£n, t), (£n+1, t), (£n+2, t) and such 
that 2,41, — £r = Az, = c is constant for any r. That is, we will restrict to lattices 
with sides of equal length. The cross-section (29) evaluated on lattices of this form 
will be given by 


1 
Pr Un =€3, P+ (Un+1 — Un)/C = €2, p: (z3 (Un+2 2un+1 + Un)) = e1 


where a = (1/c?) det(un42,Un41,Un) and u; = u(x,t). The left moving frame 
associated to this cross section is given by 
1 1 
-1 


P = (3 (Un+2 2Un41 un), č (Un41 un), Un) 
o? 0 0 
= (un42 Unyi Un) |2? t 0 
a 1 


which can clearly be gauged to 


n= (un+2, Un+1, un) . 
The multispace subspace (30), when restricted to our partially coalesced lattices, 


becomes 


(32) det(un+2, Un+1, Un) = @ 


for all n. Let us introduce one last gauge by the matrix 


1 0 0 
g= —An-1 1 0 
0 0 1 


If un+3 = k Unze + k3 Un+1 + Un, then, the discrete x-Maurer-Cartan matrix asso- 
ciated to p = ng is given by 


0 1 0 
Ka = | by. “Gye 1 
1 0 0 


where bn = k3 and an = Ko As before, the t-Maurer—Cartan matrix is given by 


—(w +78) wh rë 
Nn = Pn (On)t = uy wt ort 
n n n 
V2 We T2 
and the local basic syzygy is 
=i = 
Ky (Knj = Nn — Ka, Nn Kn. 
This syzygy solves for N,, in terms of r; 
n — „n+l. n _ „n+l n+l n+l. n — „n+l. 
wo =r; wi =r3 +anri +bnro ; w3 =TO ; 
= 2 1 
vl =r) +bn-irt; vs = rīt — anro” ; 


and it provides the condition for preserving the restriction to the multispace sub- 
manifold (32), namely 


(33) rt? + ase +r3 + anyir”? + begin t? + barot! =0. 
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If the map r3 > rgt? +r3+! +r} = (T? +7 + 1)r} is invertible (which is true if 
N + 3s for any s as shown in [46]), this condition solves for r§ in terms of rı and 
ro. The syzygy also describes (an)s and (by) to satisfy 


(an): = (+ Anbnsi)rot? — (anbara! + is) 
+ (bp + ananpi )ri t? — (bnr? + irr) + anr”? — anr} 
(bn) = ree? a (bnan ri"? tT) = (aa babna )ro t? 


— bape? + anro — (bara? + Dirk). 


If further we impose the syzygy 


1 2 

(34) (tnt + —— (Un42 — Gnting1) + Fun = 0 

bn—1 3 
then we can see that rý = ae r? =O and r} = —(T? +741) 70(7+ 1) 22 = -2 
is the solution of (33) for these choices. Then 

1 1 
(an)e ~ bn+1 is bn—2 
(bn) = a + m 


and the changes 

= cA B =~ _On+1_ 
bnbn+1bn+2 i k bnbn+1 

transform this equation in the integrable discretization of Boussinesq 


(an) = On(Bn+2 =< Bn-1) 
(Bn)t = An—1— An + Bn(Bn+1 = Brat): 


This system appears in [25]. 

It is not clear to us how to systematically connect integrable discrete systems and 
evolutions of polygons as given by (34). In the continuous case there is a general link 
between Hamiltonian evolutions at the level of the invariants and evolution of curves 
on geometric manifolds which are homogeneous of the form G/H or (G x R")/G 
with G semisimple (see [42], [43]), but the situation in the discrete case is not so 
clear in general. In particular, the syzygies (31) and (34) are not the restriction 
of the same syzygy on different points in the multispace, even when we account 
for all the different changes introduced by gauges. Still, it is widely known that 
certain evolutions of polygons result in completely integrable discrete systems (see, 
for example, [3] in the Euclidean case and [2] in more cases with further restrictions 
on the lattices). The multispace allows us to construct geometrically without the 
need to account for the limits. We include one more example along these lines and 
further use of multispace in this area will appear elsewhere. 


An = 


5.2.2. Discretization of the Sine-Gordon equation. It is well known that the Codazzi- 
Mainardi equations for Euclidean surfaces in RÌ with constant negative Gauss cur- 
vature includes a Sine-Gordon equation, a well known completely integrable sys- 
tem, that decouples from the rest of the determining equations for the surface. The 
Codazzi—Mainardi equations are simply syzygies for a well-chosen moving frame, 
hence using the multispace framework we will be able to find a discretization of the 
Sine-Gordon equation with strong geometric meaning as determining mixed lat- 
tice/surfaces with negative Gauss curvature. It is not clear to us if the discretization 
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below is completely integrable as it becomes part of a system of equations defining 
the lattice/surface, rather than decoupling to discretize Sine-Gordon individually. 
A discrete Sine-Gordon equation on lattices was also defined in [2], although the 
conditions that the authors imposed on their lattices are not impose here. Further 
study on the connection between both approaches will appear elsewhere. 

We review the continuous case first. 


5.2.3. Sine-Gordon as syzygy of Euclidean surfaces with constant negative curva- 
ture. Let G be the Euclidean group represented as the subgroup of GL(4, R) 


(35) g= E 6] 


where v € R? and © € SO(3). The group acts in R? with the standard action 
g:u = Ou +v which coincides with the one induced by g er Let u(x,y) be 


a parametrized surface and assume that x and y are normalized to measure the 
arc-length in the x and y direction. That is, assume that us| = [u| = 1. Let us 
define a moving frame through the normalizations 


g°u=0, g:Ug=e1, J+ Uy = cCosae, + sin aeg 


where a is the angle formed by uz and uy. Solving the equations we obtain that 
v = —Ou and 
OT = (us = (uy—cosau,) N) 


sin a 
where N = =+~(uz X uy) is the standard normal unit vector determined by the 
parametrization. Using the traditional notation Uss = Tł4Us +r? uy + eN, usy = 
Tius + PZouy + fN, uyy = [Que +I22Uy + gN, we can write the Maurer-Cartan 
matrices as 


oaia A O\/i 0\_ 7/0 0 
~ PrP = \-—(Ou) Or) Lu OF) ~ \-e, 0,0T 


5 0 0 
E act eee 
el = (_ cos ae; — sin aeg ee") 


where 
0 sin ar? e 
(36) N=0,07 = —sinal7, 0 sa (f — cos ae) 
e =. (f — cos ae) 0 
0 sinal 7, f 
(37) K=0,07 = |-sinar?, 0 (g — cosa f) 
-f = (g — cosa f) 0 


Substituting these values in the local syzygy K, = Ñ, + [K ; N ] and selecting the 
R component, we get the equation 


— cosa 0 0 sin ar? 
0= | -sina | — | —sinar?, | +cosa | —sinal}, | +sina 0 
=i 
0 ` —f —e ana (f — cos ae) 


The last entry is trivial, the first solves for the value 


gl ES 
r= 
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while the second one simplifies to r?, = 0, whenever sina 4 0. The SO(3) portion 
of the syzygy is given by 

(38) K, = N, + [K, N]. 

If we write down the equation that does not involve derivatives of the second fun- 
damental form (the equations defined by the (1,2) entry), we have 


(sinaTî4), = m (f — cosae) + — (g — cosa f), 
Y sma sin Q 
which becomes i 
Laut p aga 
sy e leg — f). 
If K is the Gauss curvature, we know that K = ££ 7, where E = us| = G = 


EG- F7? 
|u,|| = 1 and F = ux : uy = cosa. Thus, the equation becomes 


Qzy = —sinak 


F which is the Sine—Gordon equation whenever K is constant and negative. 

Notice that this equation is not enough to determine the surface. Indeed, solving 
for a only determines the first fundamental form (or metric), and the knowledge 
of the Gauss curvature does not suffice to determine the second fundamental form. 
Indeed, one would need two more equations to do so, given by the two remaining 
entries (1,3) and (2,3) of the SO(3) portion of the local syzygy, i.e. (38). Thus, the 
surface is determined upon solving a system of 3 equations, one of which decouples 
and is equal to Sine-Gordon. 


5.2.4. A discrete Sine-Gordon evolution as a syzygy of a mixed difference-discrete 
lattice. Assume we have a smooth family of polygons, or a mixed (1, 1)-lattice (one 
continuous direction and one discrete one) of the form yr, Yr+ı in the y direction 
and continuous 1-jet in the x direction. As far as we use the same multispace 
cross-section, we will have guaranteed that discrete invariant data approximates 
the continuous one. Thus, consider the transverse section 
(g:TUp — g: ur) = [Aur] r 
Yı — Yo Ayo 
where g is as in (35); Tur = Upy1 = u(x, Tyr) = U(X, Yr+1); *p is the prolonged 
action given by g-pte = Ouz; and wp is an unit vector with Wp = cos ane; +sin anez 
where an = Q(x) is the angle between uy and Aur, Au, = (T — 1)u,. From here, 
the multispace cross section defines v, = —O,u,, Oruz = |uz|e1 and 0,Au, = 
[Aver fw . 

With these choices the right moving frame becomes 


1 0 
Pr D i 3 , Ur = —0,Ur 


ef = (tr Np br) 
_ (Unde. cat 1 ( Au, 

"O Kurel "sina, \|Au,| 
From now on, and for convenience, we will drop the subindex to denote position 
unless the situation is confusing, indicating a change in position by the application 
of the shift operator (T*u = un). We will calculate the left Maurer-Cartan 


9° Ur =g: U(t, Yr) =0; gp Ue = |vcler; 


with 


cos arts) a bp St Xe: 
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matrices, the more geometrically significant one (those in the continuous case are 
right ones). The left Maurer—Cartan matrices are given by 


ne ee E ac e ae a A 
ESPEP =e 3) a rer) = (jaupo J 


where K = @TO7!; and by 


N = p(o) = (a g) a er,) ~ (iles N) 


where N = O(©7!)„. The local syzygies are given by 
(39) (R) =KTN — NK. 
As we did in the continuous case, and to ensure convergence, we will restrict to the 
submanifold of the multispace defined locally by 
[Aun _ 
[Ayn] 
For simplicity we will restrict further to those lattices where |Ay,| = € for all r (and 
hence |Au,|| = €). 
Equation (39) breaks into two equations, namely 


uel =1 


(40) We = Ke, —Nw-e, 
(41) K, = KTN-NK. 
Assume K factorizes as 
Q 0 
(42) K=Y a 7 
with Q € SO(2), for some Y = exp & i) (this is always possible when TO 


is closed enough to © so that K is closed enough to the identity). Assume further 


that 
cosk —sink Ls 
D= e cos k ) =) 
v and U denoting the two columns of 2. We will denote with a hat the transforma- 


tion 0 = F oe v = Jv, and so © = Jw. Finally, denote N by 


wwe) a= (4 9). (8) 


With this notation (40) can be rewritten and simplified to equations 


(44) —sy;u = zw 
(45) wey = wrvtcy-vy-w 
1 
4 p= ' 1 ; 
(46) a v aA cos k — cy: v y2) 


where s = GI sin |y|, c = pr (cos lul — 1), and as usual y? = (y1, y2). 
The remaining three equations that will determine the lattice/surface are given 
by the three entries in the s0(3) component (41). We will only reproduce the portion 
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corresponding to Sine-Gordon, that is the (2,1) entry of (41). After some long, but 
straightforward algebraic manipulations, the equation becomes 


(47) kz + cdet(y, Ye) = Tv — cos |ylv + sdet(y, z). 
Looking at (37) and (36), and comparing it to (42) and (43), we see that 
s det(y, z) 


discretizes the determinant of the first two entries of the last column in both (37) 
and (36). That is, sdet(y, z) discretizes 


f e — 
det ( 1 Gg cosas) =x (f — cos ae) = —sinak 


were K is the Gauss curvature of the surface. Therefore, we can define 


Kesta 


° = . det (y, z) 
sina lyl sina 


to be the discrete Gauss curvature for the lattice/surface. Then, (47) becomes 
ky, + cdet(y, yx) = Tv — cos |y] — sin aK. 


This will be a discretization of Sine-Gordon, together with the other equations in 
the system. 


APPENDIX A. EQUICONTINUOUS FAMILIES OF DISCRETE FRAMES 


In this Appendix, we use the Arzela-Ascoli Theorem to give a general convergence 
result for an equicontinuous family of moving frames. This provides a rigorous 
foundation to a variety of examples involving the discretisation of a smooth frame. 

Let M be a manifold, and let G be a Lie group with local metric d. The set 
G™ consists of all continuous maps from M to G, and we give it the compact-open 
topology, defined as that generated by finite intersections of the so-called subbasic 
sets, 

(A,V) ={f E€ G| f(A) CV} 


where A C M is open and V C G is compact. A sequence of maps converging in 
this topology is uniformly convergent on compact subsets. 


Definition A.1. A family F C G™ is said to be equicontinuous at yo € M if for 
all e > 0 there exists a neighbourhood U (yo) C M such that for all p € F, 


pU (yo)) C B(p(yo), €) = {9 € G| d(g, p(yo)) < €}. 


Theorem A.2. Suppose that a family of left (resp. right) moving frames F c GM 
satisfies the following: 


(1) F is equicontinuous on M, and 
(2) the set 


{ply |p © F} 
has compact closure for each y. 


Then F is a compact and equicontinuous family of moving frames. 
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Proof. We give the proof for F a family of left frames, the proof for right frames is 
analogous. The conditions of the Theorem are precisely those of the Arzela-Ascoli 
Theorem, ({15], XII, Theorem 6.4), which yields that the family F is compact and 
equicontinuous. We need only show that its elements are also equivariant with 
respect to the group action on M. Fix y € M and f € F, and let e > 0 be given. 
By the definition of the closure of the set F and the continuity of the group action, 
there is a p € F and a neighbourhood U of the identity e € G such that for g € U 


we have both d(f(g-y),e(g-y)) < 3€ and d(g- f(y), g: p(y)) < 36. Then 


d(g- f(y), f(g-y)) d(f(g-y),e(g-y)) + d(e(g-y), g: f(y)) 
= df(g-y),elg-y))+d(g- ely), g: f(y) 
< €, 


A 


so that f is equivariant, as required. 


Example A.3. Consider the scaling and translation action of R? given on a equi- 
variant family of Lipschitz continuous curves (x, y(x)) in the plane by g- (x, y(x)) = 
(x,exp(A)y + k). A smooth frame is given by g-y = 0, g-yx = 1, or (exp(A), k)I = 
(1/y2,—y/Yx). Suppose now we wish to discretise this frame in a way that is com- 
patible with the smooth frame and with forward difference, that is Yn+1 = Yn + ayz. 
Then the discrete frame pn would be obtained by the normalisation equations, 
9° Yn =Q, 9° Ynti = a, so that 


( a —AYn ) 
Pn = , : 
Yn+1 — Un Yn+1 — Yn 


This family of frames is straightforwardly seen to be equicontinuous, to have the 
smooth frame as its continuum limit, and to have the smooth Maurer Cartan in- 
variants as the limit of the discrete ones, provided the parameter a scales as the 
mesh size n41 — Ln. 
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